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ABSTRACT 

We present extensive optical {UBVRI, g'r'i'z', and open CCD) and near-infrared {ZY JH) photom¬ 
etry for the very nearby Type HP SN 2013ej extending from -|-1 to -1-461 days after shock breakout, 
estimated to be MJD 56496.9 ± 0.3. Substantial time series ultraviolet and optical spectroscopy 
obtained from -|-8 to -1-135 days are also presented. Considering well-observed SNe HP from the lit¬ 
erature, we derive UBVRIJHK bolometric calibrations from UBVRI and unfiltered measurements 
that potentially reach 2% precision with a B — V color-dependent correction. We observe moderately 
strong Si H A6355 as early as -1-8 days. The photospheric velocity (uph) is determined by modeling 
the spectra in the vicinity of Fe H A5169 whenever observed, and interpolating at photometric epochs 
based on a semianalytic method. This gives Uph = 4500 ± 500 km s“* at -1-50 days. We also observe 
spectral homogeneity of ultraviolet spectra at -1-10-12 days for SNe HP, while variations are evident 
a week after explosion. Using the expanding photosphere method, from combined analysis of SN 
2013ej and SN 2002ap, we estimate the distance to the host galaxy to be O.Oilip g Mpc, consistent with 
distance estimates from other methods. Photometric and spectroscopic analysis during the plateau 
phase, which we estimated to be 94 ± 7 days long, yields an explosion energy of 0.9 ± 0.3 x 10®* ergs, 
a final pre-explosion progenitor mass of 15.2 ± 4.2 Mq and a radius of 250 ± 70 R©. We observe a 
broken exponential proHle beyond -1-120 days, with a break point at -1-183 ± 16 days. Measurements 
beyond this break time yield a ®®Ni mass of 0.013 ± 0.001 M©. 

Subject headings: supernovae: general — supervovae: individual (SN 2013ej) — galaxies: distances 
and redshifts — techniques: — photometric 


1. INTRODUCTION 

Supernovae (SNe) exhibiting subst antial hydrogen i n 
their spectra are classified as Type H (Filippenko 1997). 
These events are considered to result from the sudden 
core collapse (CC) of massive stars that still retain sub¬ 
stantial hydrogen envelopes. Early-time spectra are ba¬ 
sically blue continua with P Cygni lines of hydrogen. 
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SNe H manifest in a variety of subtypes, with Type 
SNe HP yielding distinctive plateaus of bright optical 
emission lasting roughly 100 days. The plateau phase is 
believed to arise from a particularly extended hydrogen 
outer layer that sustains optical emission through recom¬ 
bination as the photosphere recedes and the outer enve¬ 
lope cools over time. After the plateau phase ends, sub¬ 
sequent evolution is powered by radioactive decay. This 
behavior yields direct measurement of radioactive mate¬ 
rial produced from the explosion. While some variation is 
observed in the late-time properties among SNe HP, vari¬ 
ation is more evident in the properties during early times 
and the photospheric phase, such as rise time, absolute 


peak magnitude, plateau length and slope (e.g. Ander¬ 


son et al. (2014)). Unlike thermonuclear SNe la, which 
are thought to come mostly from near-Chandrasekhar- 
mass white dwarf thermonuclear explosions, SNe HP are 


believed to arise from massive pr ogenitors (Heger et al. 


2003 


Utrobin & Chugai 2009) ranging from 8 to 25 


Using pre-SiN imaging data, Smartt et al. 


(2009) 

obtained a Zero Age Main Sequence (ZAMS) mass range 
of 8 to 17 Mq for these events. Nevertheless, their char¬ 
acteristics have lent themselves to use as cosmic distance 
indic ators and possible independent probes of dark en- 


!rgy i Hamuy et al.||2001| Hamuy fc Pinto |[2U^ Nugent 


et al. 2006 ' Poznanski et al.||2010 ). 


SiNe HP present the opportunity to measure a wealth 
of physical parameters from the explosion, and the exten¬ 
sive data available for nearby events are crucial to pin- 
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ning down the mechanisms involved. This, in turn, is im¬ 
portant to any use as cosmological probes from the most 
frequently occuring SN types (e.g., Li et al. (20111). 

On 2013 July 25 (UT dates are used throughout this 
paper), discovery with the 0.76-m Katzman Automatic 
Imaging Telescope (KAIT) at Li ck Observatory of a new 
SN IIP in M74 was announced (Kimetal. ||2013). This 
made SN 2013ej one of the closest SiNe ever discovered. 
Prediscov ery photometry was obtained with the Lulin 
telescope (|Lee et al.|20I3 ) a nd the ROTSE-IEb tel escope 
at McDonald Observatory (|Dhungana et al.|2013|), mak¬ 
ing this also one of the best-observed young SNe IIP. 
Follow-up spectroscopy was performed using the Hobby 
Eberly Telescope (HET), and the Kast spectrograph at 
Lick O bservatory, providing a classihcation and a red- 
shift. Valenti et al. (2014) performed an analysis of 
the first month ot photometry and spectroscopy, yield¬ 
ing constraints on the object and indicating it to be one 
of the more slowly evolving SNe IIP at early times. They 
identified a moderately strong Si II feature, blueward of 
Ha in the first month. Pre-explosion images obtained 
with the Hubble Sp ace Telescope (HST) were analyzed by 


Fraser et al. (2014), from which they proposed two pos- 
sible progenitors, with the redder source being the more 
likely candidate. Using an M-type supergiant bolometric 
correction, they estimated the mass of the progenitor to 
be 8-15 .5 Mf7^. More recently , from hydrodynamic simu¬ 
lations, Huang et al. (2015) found the progenitor to be 
a red supergiant with a derived mass of 12-13 Mq prior 
to explo sion. We also note that, from an independent 
dataset, Bose et al. (2015a) have favored SN 2013ej to 
be a Type HL event, accounting for the observed steep 
plateau and the systematically high velocity of strong H I 
lines. 

We present an extensive analysis of unfiltered CCD 
and broadband photometry from the ultraviolet (UV) 
through the infrared (IR), and a time series of UV and 
optical spectroscopy, for SN 2013ej. We consider all the 
measurements relative to 2013 July 23.9 (MJD 56496.9) 
unless otherwise explicitly stated. Sectionpresents the 
data, while Sectionj^ describes the photometric and spec¬ 
troscopic reductions. Utilizing open-CCD and broad¬ 
band photometry, we analyze the early-time photometry 
to derive the time of shock breakout in Section a This 
section also presents bolometric calibration of unhltered 
and broadband photometry, as well as a derivation of 
photometric observables such as color and temperature. 
Analysis of UV and optical spectroscopic features from 
+8 d to -1-135 d is discussed in Section^ In Sectionawe 
derive the photospheric velocity at photometric epochs, 
from which we utilize the expanding photosphere method 
(EPM) to estimate the distance to SN 2013ej. Kinemat¬ 
ics of the explosion, properties of the progenitor, Ni mass 
yield, and other physical properties are derived in Sec¬ 
tion The discussion and our conclusions are presented 
in Section [8l 


2. OBSERVATIONS 
2.1. Photometry 


SN 2013ej was discovered by the Lick Observatory Su¬ 
pernova Se arch (LOSS; Filip penko et al. 2001) on 2013 
July 25.45 (Kim et al. |2013), using unhltered data taken 
with KAI'i. A color combined frame of the SN 2013ej 



Fig. 1. — The field around SN 2013ej on a color-combined {BVI) 
CCD frame taken with the 0.6 m Schmidt telescope at Konkoly 
Observatory, Piszkesteto, Hungary. Konkoly photometry reference 
stars in the vicinity of SN 2013ej are shown. 


and the host galaxy M74 is shown in Fig. [H The 
0.45 m ROTSE-HIb telescope also observed SN 2013ej 
in automated sky patrol mode, hrst on 2013 July 31.36. 
ROTSE-HIb is operated with an unhltered CCD with 
broad wavelength transmission over the range 3000- 
10,600 A. Precursor ROTSE images from July 14.42 
rule out any emission at a limiting magnitude of 16.8. 
Careful analysis of additional ROTSE-HIb observations 
reveals the earliest detection at July 25.38, about 100 


minu tes prior to the discovery epoch (Dhungana et al. 


2013). Eollowing discovery, we scheduled follow-up obser- 


vations with the goal of obtaining well-sampled photom¬ 
etry of this bright, nearby SN. Unfortunately, weather 
conditions were not optimal for the following 5 days for 
ROTSE-HIb when the SN was nearing its peak. We then 
continued observations for 200 days (see Fig. [^. 

We obtained broadband photometry with th^O/90 cm 
Schmidt telescope of the Konkoly Observatory at 
Piszkesteto Mountain Station, Hungary, through Bessell 
BVRI hlters. This Konkoly dataset spans from -|-8 d to 
-1-130 d. Photometric observations were also performed 
at Baja Observatory, Hungary, with the 50 cm BART 
telescope equipped with an Apogee-Ultra CCD and Sloan 
g'r'i'z' filters. 

Photometry was also obtained with the multi-channel 
Reionization And Transients InfraRed camera (RATIR; 


Butler et al. 2012) mounted on the 1.5 m Johnson tele- 
scope at the Mexican Observatorio Astronorhico Na- 
cional o n Sierra San Pedro Martir in Baja California, 
Mexico (Watson et al. 2012). Typical observations in¬ 
clude a series ot 80 s exposures in the ri bands and 60 s 
exposures in the ZYJH bands, with dithering between 
exposures. Near-IR data from RATIR span from -|-3 d 
to -I-I25 d. 


In addition to the unhltered data taken in discovery 
mode, scheduled follow-up Bessell BVRI photometry 
was obtained with KAIT and the Nickel 1-meter tele¬ 
scope located at Lick Observatory. Starting on June 30, 
a thorough sample of both unhltered and BVRI mea- 
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TABLE 1 

Tertiary Konkoly BEJ?/ measurements of the standard stars in the vicinity of SN 2013ej used eor Konkoly photometry'^ 


ytar 

a 

(J2000) 

d 

(J2000) 

B 

(mag) 

V 

(mag) 

R 

(mag) 

I 

(mag) 

2MASS® J01365863+1547463 (A) 

2MASS J01365760+1546218 (B) 
2MASS J01365154-I-1548473 (C) 
2MASS J01364487+1549344 (D) 

01:36:58.60 

01:36:57.56 

01:36:51.51 

01:36:44.88 

-1-15:47:47.32 

-1-15:46:22.04 

-1-15:48:48.04 

+15:49:35.88 

13.19 (0.02) 
13.95 (0.02) 
14.62 (0.03) 
15.74 (0.03) 

12.60 (.01) 
13.16 (.02) 
13.99 (.02) 
14.93 (.02) 

12.32 (0.02) 
12.77 (0.02) 
13.69 (0.02) 
14.59 (0.02) 

11.90 (0.02) 
12.30 (0.02) 
13.26 (0.02) 
14.10 (0.02) 


Photometric uncertainties are given inside parentheses. 

^ Two Micron Ali-Sky Survey 

surements was obtained until late in the nebular phase. 
Unhltered KAIT data extend to +213 d, while BVRI 
data span from +7 d to +461 d (see Fig. and Fig. 

SN 2013ej was also monitored with the UVOT instru¬ 
ment onboard the NASA Swift space telescope through 
the uvw2, uvm2, uvwl, u, b, v filters. These frames were 
collected from the Swift archivj^ Swift data range from 
+ 7 d to +138 d The Swift dataset has be e n published 
bylVale nti et al.|(|2014|), [Huang et al. |(|2015|),|Bose et al. 
(2(J15a(). (Jur reduction of Swift irames is consistent with 
these works in the u, b and v hlters in the plateau, and 
until +30 d after the explosion in the uvw2, uvm2 and 
uvwl filters. Ho wever, we obtain sig nificantly brighter 


(2015) beyond +30 d in 
(diven fhe tact that uvw2, and 


magnitudes than pTuangetaL 
the later three filters. 
uvwl filters have extended red tails (e.g. Ergon et al. 
2014) and using our photometry, all three of these have 


a marginal contribution to the total flux after +30 d 
we ignore the flux from uvw2, uv m2 a nd uvwl bands 
beyond this epoch (also see Section 14.5). 

We also note that the detection of SN 20 13ej at its 
young est observed phase was announced by |Lee et ^ 
(2013) in the BVR bands, on July 24.8, which is 15 hr 
earlier than the first ROTSE-IIIb detection. A nonpho¬ 
tometric prediscovery detection on images taken on July 
24.125 was also reported by C. Feliciano on the Bright 
Supernovae websiteP^ No emission on 2 013 July 23.54 
at E = 16.7 mag was seen by ASAS-SN (Shappee et al. 

Ml. 


2.2. Optical and Ultraviolet Spectra 

A total of 17 low-resolution optical spectra of SN 
2013ej were obtained using the Marcario Low-Resolution 
Spectrograph (LRS; |Hill et al.| (|1998|)) on the 9.2 m 
Hobby-Eberly Telescope (MET) at iVlcUonald Observa- 
tory, the dual-arm Kast spectrograph (Miller & Stone 


1993 1 on the Lick 3 m Shane telescope, and t he UEep 
Imaging M ulti-Object Spectrograph (DEIMOS; Faber et 


al. (20031) on the Keck-II 10 m telescope. Fhe Kast 
and DEIMOS observations were aligned alon g the paral¬ 
lactic angle to reduce differential light losses (Filippenko 


||1982 ). These optical spectra span from +8 d to +135 d. 

JNear-UV spectra of SN 2013ej were taken with 
UVOT/UGRISM onboard Swift, covering the wave¬ 
length range 2000-5000 A and spanning +8-16 d. 

3. DATA REDUCTION 
3.1. Photometry 

http://heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/swift.pl 
http://www.rochesterastronomy.org/supernova.html 


reduction pipeline ( 

Yuan & Akerlof 

2008), followed by 

a DAUFMUi-based pomt-si 

)read-tui 

iction (PSl j pho- 
Because of signih- 

tometry technique 

Stetson 

19871). 


age differencing, we performed aper ture photometry of 
SN 2013ej (e.g., |Marion et al.| (|2015|)). An aperture size 
of 1 full width at half-maximuin intensity (FWHM) of the 
median PSF on each image was considered, and we chose 
a background-sky annulus having inner and outer radii 
of 2 and 4.5 times the FWHM. Additionally, a reference 
template image was smeared to reflect the PSF at each 
epoch, and the underlying host-galaxy contribution in¬ 
side the aperture was subtracted. The typical FWHM of 
the PSF during the observation timescale was 3-4". We 
calibrated the derived relative flux to the R band from 
the USNO Bl.O catalog. The instrumental calibration 
and comparison to other data for analysis are presented 
in Section ID 

Filtered data from Konkoly were reduced with stan¬ 
dard IRAlf^ routines to get the SN magnitudes. The 
instruments magnitudes were transformed to the stan¬ 
dard J ohnson-Cousins system via local tertiary standards 
tied to Landolt (1992) standards on a photometric night 
(see Table 'jTpina "Fig. 1^. The g’r’i’z’ data from Baja 
Observatory were standardized using ~ 100 stars within 
the ~ 40 X 40 arcmin^ field of view around the SN, taken 
from the Sloan Digital Sky Survey (SDSS) Data Release 
12 catalog. In order to avoid selecting saturated stars 
from the SDSS catalog, a magnitude cut 14 < r' < 18 
was applied during the photometric calibration. 

PSF photometry was performed on KAIT and 
Nickel reduced data (|Ganeshalingam et al.||2010| using 
DAOPHOT. Several nearby stars were chosen from the 
APAS^^ catalog, and the rnagnitudes were first trans¬ 
formed to the Landolt systenp^ before calibrating KAIT 
data. We used APASS R band magnitudes to calibrate 
the KAIT unhltered photometry. Image subtraction was 
not performed for KAIT data, as the object was ex¬ 
tremely bright and far from the galaxy core. 

For RATIR data reduction, no off-target sky frames 
were obtained on the optical CCDs, but the small galaxy 
size and sufficient dithering allowed for a sky frame to be 
created from a median stack of all the images in each 
hlter. Flat-held frames consist of evening sky exposures. 


IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation (NSF). 
http://www.aavso.org/apass 

http://www.sdss.org/dr7/algorithms/sdssUBVRITransform.html 
+Lupton2005 
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KAIT I-1.5 
Konkoly 1-1.5 
KAIT R-1 
Konkoly R-1 
KAIT Unf. 
ROTSE Unf. 
KAIT V-i-0.5 
Konkoly V-i-0.5 
KAIT B-i-1 
Konkoly B-i-1 
Swift v-i-3.0 
Swift b-i-3.0 
Swift u-i-3.0 
Swift uvwl -1-3.0 
Swift uvm2-i-3.0 
Swift UVW2-I-3.0 
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100 150 

MJD-56496.9 



Fig. 2.— Open CCD and multiband photometry of SN 2013ej. FiG. 3.— SDSS g’r’i’z’ photometry from Baja Observatory, as 

KAIT BVRI and unfiltered data points are shown with empty well as RATIR optical and near-IR photometry, of SN 2013ej. 

circles, while Konkoly BVRI and ROTSE points are solid circles. 

Swift photometry is represented with filled square symbols. 


Given the lack of a cold shutter in RATIR’s design, IR 
dark frames are not available. Laboratory testing, how¬ 
ever, confirms that the dark cu rrent is negligible in both 
IR detectors (Fox et al. 2012). RATIR data were re¬ 
duced, coadded, and analyzed using standard CCD and 
IR processing techniques in IDL and Python, utilizing 
online astrometry programs SExtractor and SWarpp^ 
Calibration was perform ed using field stars with reportec 
fluxes in both 2MASS (Skrutski e et al. 2006|) an d the 
SDSS Data Release 9 Catalogue (Ahn et al. 2012). 


Figure shows the final calibrated SiN zUIiJej light 
curves in the ROTSE and KAIT unfiltered bands, the 
KAIT and Konkoly BVRI bands, and the Swift UVOT 
bands. Comparison of the data from various sources re¬ 
vealed that they are generally consistent within ±0.1 
mag in all optical bands. Figure [^illustrates g’r’i’z’ 
Baja photometry and riZY JH RATIR photometry. The 
Appendix tabulates ROTSE, Konkoly, Baja, KAIT and 
Nickel, and RATIR photometry in Tables ll\ 12 13 14 
and[T^ respectively. 


TABLE 2 

Observing log of SN 2013ej optical spectra. 


UT Date 

MJD Epoch'^ (days) Instrument 

2013 Aug 1.41 

56505.41 

-t8 

HET 

2013 Aug 2.46 

56506.46 

+9 

DEIMOS 

2013 Aug 4.38 

56508.38 

-til 

HET 

2013 Aug 4.51 

56508.51 

-til 

Kast 

2013 Aug 8.52 

56512.52 

-tl5 

Kast 

2013 Aug 12.50 

56516.50 

-tl9 

Kast 

2013 Aug 30.50 

56534.50 

-t37 

Kast 

2013 Sept 6.41 

56541.41 

4-44 

DEIMOS 

2013 Sept 10.60 

56545.60 

4-48 

Kast 

2013 Oct 1.54 

56566.54 

4-69 

Kast 

2013 Oct 5.33 

56570.33 

4-73 

Kast 

2013 Oct 8.48 

56573.48 

4-76 

DEIMOS 

2013 Oct 10.29 

56575.29 

4-78 

Kast 

2013 Oct 26.26 

56591.26 

4-94 

Kast 

2013 Nov 2.34 

56598.34 

4-101 

Kast 

2013 Nov 8.31 

56604.31 

4-107 

Kast 

2013 Nov 28.37 

56624.37 

+127 

Kast 

2013 Dec 6.39 

56632.39 

+135 

Kast 


^ Epochs are rounded to days since explosion. 


3.2. Spectroscopy 

All of our opti cal spectra were reduce d using standard 
techniques (e.g., Silverman et al. (2012)). Routine CCD 


processing and spectrum extraction were completed with 
IRAF, and the data were extracted with the optimal al¬ 
gorithm of 


Horne 


(1986). We obtained the wavelength 
scale from low-order polynomial fits to calibration-lamp 
spectra. Small wavelength shifts were then applied to the 
data after cross-correlating a template sky to the night- 
sky lines that were extracted with the SN. Using our 
own IDL routines, we fit spectrophotometric standard- 


star spectra to the data in order to flux calibrate our 
spectra and to remove tellur ic lines (Wade & Horne 


1988 Matheson et al. 2000). A log of observed opti¬ 


cal spectra is given m Table and plotte d in Figure 


HET spect ra are archived on WISEREI 
Yam||2012 ), and all of our spectra will 


22 


(Yaron & Gff 
^e made publicly 

available from the database. 

UV spectra were collected from the Swift archive, 
and were reduced using the uvotimgrism task in 
HEAsof^^ The log of the UGRISM spectral obser¬ 
vations is given in Table and the spectra are plotted in 


SExtractor and SWarp can be accessed from 
htt p: / / WWW. astromat ic. net / software. 


http://wiserep.weizmann.ac.il 
http://heasarc.nasa.gov/lheasoft/ 
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TABLE 3 

Observing log of Swift UVOT/UGRISM spegtra 


UT Date 

MJD 

Epoch 

(days) 

Exposure 

(s) 

S/N 

2013 Jul 31.8 

56504.8 

-t8 

4142 

30 

2013 Aug 3.1 

56507.1 

-tio 

4946 

41 

2013 Aug 4.8 

56508.8 

-tl2 

4896 

36 

2013 Aug 7.4 

56511.4 

-tl4 

3861 

25 

2013 Aug 8.2 

56512.2 

-tl5 

4449 

21 

2013 Aug 9.2 

56513.2 

-tl6 

4449 

25 



Rest Wavelength(A) 

Fig. 4.— Time series of optical spectra of SN 2013ej. Phases are 
in days since e™losion (MJD 56496.9). The log of observations is 
given in Table 1^ 


Figure 


4. PHOTOMETRIC ANALYSIS 
From the lack of narrow Na I D lines from the host 


galaxy, Valenti et al. (2014) showed that the reddening 
from M74 in the direction toward SN 2013ej is negligible. 
No evidence of Na I D lines from the host was seen in 


spectra of 
Thus, we c 
the Milky 

Bose et al. ( 

20158 

l) and in our own sample. 

Lo not consider any host extinction. We adopt 
Wav reddening value o^ E(B — V) = 0.061 mag 

(Schlafly & Einkbeiner 

2011 

in our data sample. We 
used E{B — V)tot = 0.12 

note that Huang et af. 

2015 


while Bose et al. (2015a) adopted E{B — V)tot = 0.06 
mag. 


4.1. Early-Time Photometry 

After explosion, the gravitational waves and neutrinos 
soon escape while the electromagnetic signal is initially 
trapped in the envelope. Only when the hydrodynamic 
front reaches the photosphere, which takes hours to days, 
is the rise in intensity from the star observed. This epoch 
indicates the time of first light and marks the beginning 
of the shock-breakout phase. Precise knowledge of the 



Rest wavelength (A) 

Fig. 5.— The spectral evolution of SN 2013ej in the UV based 
on Swift/UGKISM observations. Labels next to each spectrum 
indicate the days sin ce explosion fMJD 5 6496.9). Feature identifi¬ 
cations are based on |Bufano et al^ | |2009| l. 


epoch of shock breakout is crucial to constrain explo¬ 
sion parameters and progenitor properties. It is also in¬ 
strumental for distance estimates using methods such as 
EPM (see Section]^. 

Several efforts have been carried out to model the shock 
breakout of th e compact progenitor o f SN 1987A (e.g.. 


Hdflich (1991); Eastman et al. (1994)). Notably, it has 
been shown that the breakout pe ak depends upon enve - 
lope mass and density structure ( Ealk fc Arnett p977 ), 
so the very early light curve may provide clues on the en- 
velope structure of massive stars. Recently, substantial 
theoretical studies have been performed with the goal of 
understanding the shock breakout of SNe II through a 
variety of proc e sses in several progenitor scena rio s (e.g., 
Nakar & Sari ( 2010| ; Svirski fc Nakar 1 (2014)). Couch 


et al. 1 2U15) argue that strong aspherical shocks can lead 


to breakout at different times along the periphery of the 
star compared to a spherical shock from a spherical star. 

To estimate the time of shock breakout for SN 2013ej, 
we consider several datasets during the first few days. To 
study the rise behavior, we combine ROTSE and KAIT 
unfiltered data, calibrated to R magnitudes, with the ear¬ 
liest prediscovery R band detection from Lulin Observa¬ 
tory. Because of the lack of sufficient data points in any of 
the independent datasets, we calibrate the Lulin R mag¬ 
nitude to the ROTSE magnitude in the following way. 
We saw a systematic variation of KAIT and Konkoly R 
band photometry of SN 2013ej. Allowing a similar offset 
to exist between Lulin R and KAIT or Konkoly R, we cal¬ 
culate the average of differences of KAIT R and Konkoly 
R magnitudes with ROTSE unfiltered magnitudes and 
add this as a correction to the Lulin data point to bring 
it to the ROTSE system. Eor this, we limit the observa¬ 
tions to between -1-30 d and -1-90 d in the plateau, where 
they are densely sampled in both sets and the spectral en¬ 
ergy distribution (SED) is smooth compar ed t o the rapid 
evolution during early times (see Section [4?^ . Eurther- 
more, we add an additional systematic uncertainty to the 
Lulin point from the root-mean square (RMS) of KAIT 
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R and Konkoly R magnitude differences. We note that 
the Lulin observation already has an uncertainty higher 
than 0.2 mag. Any systematic uncertainty that arises 
from translating the calibration from plateau to early rise 


time is likely to be smaller than this. Specifically, But¬ 


ler et al. (2006) found a correction of around of 0.1 mag 
between KAl'i' unfiltered and Lulin R band magnitudes 
for Gamma Ray Burst study. ROTSE unfiltered and 
KAIT unfiltered data, both of which closely track the 
R magnitudes to early times, have been independently 


time studies of both SN la and SN IIP (e.g. 

Quimby et 

al. 12007 

Zheng et al.|20I3 

). On the first nig 

tit of detec- 


about 2 hr between coadds of two sets of images. As the 
SN was still young (~ 1 d after explosion), a significant 
rise even in only 2 hr is detectable. 

A functional form of the SN IIP rise behavior has not 
been well established. A simple power law, specifically 
a rise law, has been tested in the context of SNe la 
many times, w hile recent studies (e.g. Zheng et al. 12013" 


Marion et al. 2015|) have shown departure from t^rise 
at very early times. In SNe la, the heat loss due to 
cooling of the ejecta can be thought to be compensated 
by the radioactive heating, thus maintaining the steady 
temperature, while in SNe IIP, the adiabatic cooling of 
the shock-heated envel ope is expected t o resul t in a steep 
drop of temperature. Quimby et al. (2007) fitted the 
early rise of SN IIP 2006bp with a law at very early 
times. While may be a valid approximation until a 
few days after explosion in SNe IIP, it is clearly not valid 
as long as it seems to hold in SNe la. 

Keeping this uncertainty in mind, we perform a least- 
squares fit of the rising light curve of SN 2013ej to a 
single power law, given by 


F{t) = A {t-tof, 


( 1 ) 


where A is a constant, to is the time of shock break¬ 
out, and /? is the power-law index. Only data points 
earlier than -|-2 d since explosion are considered for fit¬ 
ting, effectively including the first 4 points. None of the 
fits including data beyond -|-2 d were consistent with the 
detection and nondetections. We perform the to estima¬ 
tion relative to the Lulin observation point on July 24.8 
(MJD 56497.8). As SN 2013ej is a SN IIP with partic¬ 
ularly early photometry, we first test the power-law hy¬ 
pothesis by letting /3 float. This yields to = —2.19 days 
and /3 = 4.83 with x^/dof = 2.80 (see Fig. |^. Keeping 
/3 = 2 fixed, we obtain to = —0.90 ± 0.25 aays, which 
corresponds to July 23.9 ± 0.25. The reasonable x^/dof 
value of 1.47 indicates consistency of the early evolution 
with the t^ model until ^ 2 d after explosion. Devia¬ 
tion of /? from 2 might be indicative of asymmetry in the 
explosion itself and is still an open question. This is an 
important question for SNe IIP and requires further ob¬ 
servation of very early times. For SN 20I3ej, while the 
sparse data do not rule out the the power index of 4.83, 
the model yields better x^ /dof and is consistent with 
all reported early detections and nondetections. Noting 
this, we take MJD 56496.9 ± 0.3 as the epoch of shock 
breakout. 


4.2. Unfiltered and Broadband Photometry 



Fig. 6.— Early rise behavior of SN 2013ej. Multiple datasets 
calibrated to ROTSE magnitudes (see text) are shown. Solid lines 
are power-law fits to data obtained before +2 d since explosion. 
The triangle point is a nondetection limit on July 23.54 at V ^ 
16.7 mag, shown here for reference. The dashed line indicates the 
detection on July 24.125 with no photometry available. The inset 
illustrates the projection of where the emission would be for floating 
index (green) and fixed index 0 = 2 (blue). 


SN IIP light curves have a unique signature. After the 
shock breakout, the hot ejecta are believed to expand 
violently. The photon diffusion timescale being much 
longer than the expansion timescale, very little photon 
energy gets diffused. The ejecta would follow a homolo¬ 
gous adiabatic expansion, cooling quickly from the out¬ 
side. Soon after the ejecta cool to ^ 6000 K, hydrogen 
ions start to recombine, the opacity plummets, and dif¬ 
fusion cooling becomes dominant. This will result in an 
ionization front that recedes inward as a recombination 
wave, giving a characterstic, slowly declining, almost lin¬ 
ear, plateau phase that lasts for approximately 100 days. 
This plateau is observed as a result of decreasing opac¬ 
ity because of less scattering due to declining electron 
density. The photosphere remains contiguous with the 
receding ionisation front. After all the hydrogen recom¬ 
bines, the photosphere recedes into the inner, heavy ele¬ 
ment core and the light curve transitions to the radioac¬ 
tive tail phase. This tail is expected to decline by the 
®®Co —Fe decay at a rate of 0.98 mag per 100 days if 
all the gamma rays and positrons are trapped. 

Figure ^ shows the apparent magnitude light curves 
of SN 2013ej with unfiltered and BVRI broadband ob¬ 
servations. Each BVRI set consists of data that starts 
from around the peak, extends through a characteristic 
plateau phase lasting about 100 days, and proceeds to the 
well-observed radioactive tail phase. In the ROTSE light 
curve, the peak for SN 2013ej occurs at -1-18 d, where the 
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absolute magnitude reaches —17.5. This peak is consis¬ 
tent with the KAIT unfiltered data, both in phase and 
magnitude. On both the KAIT and Konkoly BVRI light 
curves, the peak occurs on -1-12.5 d in B, -1-15.5 d in V, 
-1-19.5 d in R, and -1-20 d in /. We do no t observe any 
obviou s secondary peak like that seen by |Bose et al. | 
(2013) in SN 2012aw at about -1-50 d in the V band, or 
an obvious minimum around -1-42 d in V, which would be 
indicative of the end of free adiabatic cooling. It is thus 
more challenging to ascertain the advent of the plateau 
phase in the photometry. We will estimate the plateau 
length in Sectionj^ From their respective peaks, the light 
curves decline by 0.038, 0.021, 0.016, and 0.012 mag per 
day in B, V, R, and I (respectively) until -1-90 d. The 
Konkoly and KAIT data are consistent in this decline 
behavior. Our V band slop e is steeper compared to the 


0.017 mag day ^ given by Bose et al. (2015a). This 


could possibly be a sampling issue, because their pho¬ 
tometry during the plateau is sparsely sampled and the 
peak is less well constrained. 

SN 2013ej has one of the steepest plateaus among 
SNe IIP (see Fig. for a comparison with other nor¬ 
mal SNe IIP). We note that the B band decline for SN 
2012aw was 1.74 mag until -1-104 d, while the R band 


showed no c hange in brightness over the plateau (Bose 
et al. |[2013). Classic SN IIP 1999em also evolved sim- 


ilarly (]Leonard et al. 2002), while the more energetic 
SN 2004et had a taster d ecline of ^ 2.2 ma g from the B 
band peak until -1-100 d (|Bose et al. |2013 ). The decline 
rate for SN 2013ej is consistently higher in all bands. 
Example SN IIL lig ht curves of the recent SN 20 13by 
( Valenti et al.||20 15) and the archetype SN 1980K (Bar- 


bon et al. 1982) are also shown. The magnitude fall of 


SiN 2013ej in the V band from peak to -1-50 d is about 
0.75 m ag. This puts SN 20 13ei within the SN IIL cate¬ 
gory of|Faran et al. | (|201 4|), where they use a cut of 0.5 
mag for a SiN IIL event. Valenti et al. (2015) observed 
a fall of 1.46 ± 0.06 mag in V tor SN 20 13by, and also 
pointed out the SN IlP-like behavior in its light-curve 
drop to the radioactive phase. They show a handful of 
objects for which the difference of the V band peak and 
-1-50 d magnitude is more than 0.5 mag, w hich would be 


in the SN IIL class of Faran et al. (2014). 


SN 2013ej, in spite of having such a steep plateau, 
also exhibits a drop at the end of the plateau that is 
significantly sharper than the decline in the plateau, 
which is also a characterstic feature of SNe IIP. A steep 
plateau for a SN IIP object may also indicate an inef¬ 
ficient thermalization of the ejecta. Additionally, with 
such a ste ep plateau, very little nickel yield might be 


expected. Bersten et al. (2011) showed from hydro- 
dynamic modelling that extensive mixing from ^®Ni is 
required to reproduce flat plateaus. The higher the Ni 
yield is, the sooner the plateau starts to flatten, and this 
will also affect the extension of the plateau duration be¬ 
cause of radioactive heating. 

When the plateau ends at about 100 d after the ex¬ 
plosion, the light curve suddenly transitions into the ra¬ 
dioactive tail. From the luminosity derived from radiac- 
tive decay of synthesized materials, in Section we will 
estimate the mass of nickel produced. We represent the 
decline behavior by separate linear fits to the data from 
-1-120 d to -1-183 d and from -1-183 d to -1-461 d. The 


TABLE 4 

Radioactive Tail Decline Behavior, x^/dof ah-E given inside 
PARANTHESIS. 


Band 

+120 d to +183 d 
(mag/100 days) 

+183 d to +461 d 
(mag/100 days) 

KAIT B 

KAIT V 

KAIT n 

KAIT I 

KAIT unfiltered 
ROTSE unfiltered 
KAIT BVRI 

1.15 ±0.08 (0.37) 
1.46 ±0.04 (1.01) 
1.54 ± 0.04 (2.66) 

1.63 ±0.03 (0.78) 
1.51 ± 0.02 (0.88) 

1.64 ± 0.07 (1.13) 
1.36 ±0.09 (0.48) 

0.75 ±0.02 (1.75) 
1.10 ±0.02 (2.69) 
1.30 ±0.01 (3.92) 
1.18 ±0.02 (2.96) 

1.06 ±0.02 (2.15) 



Fig. 7. — Comparison of some SN IIP and SN IIL light curves 
in the V band, except SN 2006bp (ROTSE). SN 2013ej has a sys¬ 
tematically steeper plateau, with a sharp drop at the end of the 
plateau. SNe IIL 2013by and 1980K have steep linear evolution 
after peak, but the drop at the end of the plateau is not as sharp. 
SN 2013ej also exhibits a steeper tail decline than the events shown 
here. 


epoch -1-120 d was chosen to ensure the late-time decay 
phase, and -1-183 was chose n as the break time of the late 
time behavior (see Section 7.1). Table lists the decay 
rate along with per degree of freedom of the respec¬ 
tive fits. It is clear that the light-curve decline in the 
tail is much steeper in all bands and unfiltered photom¬ 
etry before -1-183 d. Only B band has a slope shallower 
than ^®Co —)■ ®®Fe after -|-183d. While SN 2006 bp had a 
tail phase dec line of 0.73 ± 0.04 mag per 100 d (Quimby 


et al. 2007), which is less steep than full trapping of 


gamma rays from radioactive decay, SN 2013ej exhibits 
the opposite behavior. 


4.3. Color and SED Evolution 

The color evolution of SN 2013ej in the optical exhibits 
a rapid change in the first 30 days, as shown in Figure 
This is due to the fact that the U and B fluxes decline 
rapidly at early phases. While the evolution oi B — V 
is more rapid in the first 30 days, the V — R and V — I 
colors are smooth and slowly rising. Soon after, when the 
tem perature has fallen to around 6000 K (see Section 
4.4), the trends are more alike, and the SED is more 
uniform. Later, as the light curve approaches the tail, 
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MJD-56496.9 

Fig. 8.— Optical and near-IR color evolution of SN 2013ej. Here, 
the J — H and Y — Z colors from RATIR are shown in the Vega 
system for consistency. 


both the V — R and V — I colors show a rapid rise, as an 
effect from a greater decline of flux in V relative to the I 
band. The transition from plateau to the tail is evident 
in both optical and near-IR colors. 

In Figure]^ the evolution of the SED {\F\) is shown, 
together witn some of the contemporaneous UV and op¬ 
tical spectra (see Section!^. This observed SED evolu¬ 
tion is in agreement with rhe general characteristics of 
SNe IIP: a strong decline of the UV flux accompanied 
by a monotonic decrease of the continuum slope in the 
optical during the plateau phase, in accord with the con¬ 
tinuously reddening optical colors seen in Eigure|^ 


4.4. Photospheric Temperature 

We determine the temporal evolution of the photo- 
spheric temperature by fitting the KAIT and Konkoly 
BVI fluxes [R fluxes are omitted to avoid contamination 
from the strong Ha feature) to a Planck function at each 
epoch until ^-l-IOO d. Beyond this, the light curve enters 
the radioactive phase, and the energy mostly comes out 
in strong nebular lines. The BVI set covers the wave¬ 
length range 3935-8750 A. No UV flux is considered, as 
this would heavily bias the blackbody fits because of the 
many metallic blends occurring at shorter wavelengths. 
The temperature drops from 12,500 K at -1-8 d to 6400 K 
at -1-24 d, and it declines very slowly to 4000 K at -f 100 d 
as shown in Fig ure |10[ An independ ent estimate of the 
temperature by Valenti et al. (2014) is also shown, and 
it exhibits reasonable agreenient with our result. The 



Rest wavelength (micron) 


Fig. 9. — The evolution of the spectral energy distribution (SED) 
of SN 2013ej. Fluxes from broadband photometry from the UV, 
optical, and near-IR are plotted with filled circles, and the horizon¬ 
tal bars represent the FWHM of each filter. Phases are coded by 
colors and indicated in the legends. The optical and UV spectra at 
certain epochs (where available) are also overplotted for compari¬ 
son. Dereddening with E[B — V) = 0.061 mag has been applied. 



Fig. 10.— Evolution of the photospheric temperature using the 
KAIT (green points) and Konkoly (blue points) BVI datasets. 
Temperature estimates from [Valenti et al.|(|2014|| are shown with 
red points. 


rapid temperature drop in the first few weeks encapsu¬ 
lates quick adiabatic cooling, while later in the plateau 
phase the smooth slow decline signifies the cooling from 
photon energy diffusion during recombination at nearly 
constant temperature dictated by atomic physics. 

4.5. Bolometry 

Bolometric photometry permits determination of sev¬ 
eral explosion parameters, including the direct estimate 
of the amount of Ni synthesized during the explosion. 
UV flux in SNe la and SNe Ib/c is a small fraction of 
the total flux, since the high-energy photons are nearly 
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TABLE 5 

B — V Dependent Pseudo-Bolometric BVRI and UBVRI 
Calibration of ROTSE and KAIT Unfiltered Data 


Calibration 

Intercept 

Slope 

X^/dof 

ROTSE Unf. 

- Konkoly BVRI 

0.362 ± 0.004 

0.1004 ± 0.004 

0.54 

KAIT Unf. 

- KAIT BVRI 

0.401 ± 0.007 

0.082 ±0.006 

1.58 

RU'l'SE Unt. 

- UBVRI 

0.299 ±0.009 

0.161 ± 0.007 

1.05 

KAIT Unf. 

- UBVRI 

0.363 ± 0.006 

0.112 ±0.006 

2.71 


entirely absorbed by transition lines of ionized heavy el¬ 
ements. In SNe II, however, the UV flux dominates at 
early times. Bright UV and X-ray emission flashes are 
expected as the shock breaks out, followed by a post¬ 
break UV plateau lasting a few days. The spectral index 
evolves rapidly in the first few weeks, with a large por¬ 
tion of the flux yield in UV. Later evolution transitions 
to dominant emission in the R, I, and near-IR bands. 
Without data in all wavebands, it is generally impossible 
to obtain an exact bolometric flux. The availability of ex¬ 
tensive sets of data spanning from UV to IR wavelengths 
for SN 2013ej gives an ideal opportunity to obtain the 
most accurate estimate of the bolometric flux for SN IIP. 
Additionally, this helps to derive bolometric calibration 
relations for a broadband sample limited in wavelength 
and unfiltered sets such as from ROTSE or KAIT. We 
adopt the distance to the SN to be 9.0to;6 (see Sec¬ 
tion 1^ to estimate the integrated luminosity. Pseudo- 
bolometric and bolometric light curves of SN 2013ej are 
shown in the top panel of Figure 11 To estimate the 
bolometric flux in the late time when we do not have UV 
and NIR observations, we multiply the late time BVRI 
flux by a scale factor, which we derive by dividing the 
bolometric flux by the BVRI flux between -1-120 d and 
-1-137 d. We also note that the flux from uvw2, uvm2 
and uvwl bands contribute a total of 1% or less to the 
bolometric flux beyond -|-30d according to our red uction 
(see Fig. 11). Photometry by Huang et al. (2015) would 
contribute much less. Given the marginal contribution 
from these three bands after -1-30 d and pote ntial com¬ 


plication of red leaks for uvw2 and uvwl (e.g. Ergon et 
^|2014|, these three filters were omitted beyond -1-30 d 


m the bolometric flux calculation. 

As a first step, from the fact that the open CCD 
transmission is broad, we establish a calibration rela¬ 
tion for the ROTSE and KAIT unfiltered flux with in¬ 
tegrated BVRI flux as follows. Both BVRI datasets 
are converted to abso l ute flu x using the relations given 
(1998) corresponding to an AO 
The value of 


by Bessell et al. 
star (see 'I'able A4 of their paper). 


log]^Q(LirOTSE/ABy_R/) shows a direct relation with the 
B — V color. A linear fit of logig(LROTSE/TByB7) versus 
B — V \s shown in the top panel of Figure [T^ From the 
RMS of the residuals, we obtain a calibration precision 
of better than 5%, while about 8% precision is obtained 
from the residual without accounting for the B — V de¬ 
pendence. A similar analysis for KAIT unfiltered and 
KAIT BVRI datasets yields 4% and 6% precision, re¬ 
spectively, as shown in the bottom panel of Figure [I^ A 
summary of the fits is given in Table 
In the second step, since ROTSE aim KAIT unfiltered 




Fig. 11.— Top Panel. Integrated BVRI., UBVRI, 
UBVRIJHK, and bolometric light curves of SN 2013ej. UBVRI 
and UBVRIJHK shown here are derived using Swift u, KAIT 
BVRI, RATIR JH, and estimated flux from the K band. The 
bolometric light curve includes an additional contribution from UV 
flux below 3200 A before +30d. Beyond +30d, UBVRIJHK 
closely resembles the bolometric light curve in the plateau. Bottom 
Panel: UV, optical, and near-IR fractional flux. The UV covers 
flux below the U band, the near-IR (NIR) covers flux above the I 
band, while optical covers in between (UBVRI). The UV fraction 
is shown curtailed at the calculated end of plateau, because of the 
marginal contribution and possible contamination from red tails 
(see text). 


are also sensitive to the near-UV, we look at the be¬ 
havior by integrating UV data from Swift. As pseudo- 
bolometric flux based on Johnson-Cousins UBVRI fil¬ 
ters is commonly derived, we first ca librate Swift u to 
the Johnson-Cousins U band following Poole et al. (20081 
and integrate with the BVRI dataset. We limit the inte¬ 
gration to the wavelength range 3285-8750 A, where we 
have extended the lower and upper bounds by the half 
width at half-maximum intensity (HWHM) in the U and 
I bands. As before, from observations of the evolution of 
VglQ{L^^OTSE/L ubvri) , the RMS of the residuals after 
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Fig. 12.— Top Panel: Pseudo-bolometric BVRI calibration 
of SN 2013ej from ROTSE unfiltered photometry compared to 
Konkoly BVRI photometry. Residual from a. B — V color- 
dependent correction (histogram in blue) shows less than 5% un¬ 
certainty. The histogram shown on the left (green) is obtained 
from the residuals by comparing the two fluxes without any color 
correction. Bottom panel: Same as in top panel, but for KAIT 
unfiltered to KAIT BVRI calibration. We obtain 6% residuals by 
direct comparison and 4% residuals when applying a, B — V de¬ 
pendence. Fit equations for B — V dependence are given in Table 

El 


subtracting the mean reveals 13% precision. Calibration 
with B — V dependence improves the precision to 6% as 
shown in Figure [l^ where a x^/dof of 1.05 is obtained 
for the fit. Analogously, KAIT unfiltered to UBVRI 
(using Swift u and KAIT BVRI) yields about 5% pre¬ 
cision after a color-dependent correction, but with larger 
X^/dof. 

It is very unusual to obtain a consistently complete 
set of data for a single object in all bands and still have 
minimal systematic effects. Various calibration and cor¬ 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

B-V mag 

Fig. 13.— Pseudo-bolometric UBVRI calibration of SN 2013ej 
from the ROTSE luminosity. The B—V color-dependent correction 
(histogram in blue) improves the measurement by about a factor 
of 2 compared to direct comparison (histogram in green). 


rection methods have been developed to better estimate 
the bolometric flux, but all are limited in one way or an¬ 
other. Here we revisit this problem based on the most 
extensive sets of data from the literature. The calibra¬ 
tion sample is provided in Table 1^ This set includes 
extensive photometry from the U v to the IR. To avoid 
any confusion, we dub the values obtained by integrating 
fluxes from data as ^‘UBVRIJHKf while we label an 
equivalent flux derived from our calibration as “UtoK.” 
The same convention also holds in other cases. By “bolo¬ 
metric” flux, we mean integration from Swift uvw2 at 
blue end to H band in the red end, added with contribu¬ 
tion from K band as estimated below. We note that the 
IR flux past K band will be significant as the SN cools 
over time. We have not accounted for any correction 
from beyond K band in this procedure. 

The light curves shown in the top panel of Figure 
E] are derived by integrating data in the UBVRIJHK 
bands, the wavelength range 3285-23,850 A. To obtain 
the UBVRIJHK flux of SN 2013ej, we integrate the ob¬ 
served flux in the u band from Swift (after calibrating to 
Bessell U), BVRI from the Konkoly or KAIT data, and 
near-IR JH data from RATIR, where they are linearly 
extrapolated in the tail phase. We add an additional 
contribution from the K band by estimating the average 
fractional flux in K with respect to the UBVRIJH flux 
using the calibration sample given in Table We find 
that the K ban d contributes ^ 2% at -1-10 a, rising to 
5-6% at +80 d |Huang et al. | (|2015|). have published 
K band data but their data is rather sparse. Compar¬ 
ing their K band measurmemt with our estimated flux 
at matching epochs yielded an offset of less than 1% of 
the total bolometric flux in the plateau while they both 
agreed in the tail. 

From Figure [Tl] it is clear that the pseudo-bolometric 
UBVRI flux is significantly lower compared to the bolo¬ 
metric flux, and the difference monotonically grows over 
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TABLE 6 

SN IIP Bolometric Calibration Sample based on Well-Sampled Photometry from U through K 


Object 

Host 

Distance (Mpc) 

Total E{B — V) (mag) 

V Plateau Slope (mag/100 days) 

Feature 

References 

SN 1999em 

NGC 1637 

11.7± 1.0 

0.10 

0.31 ± 0.05 

Normal 

1,2,3,4 

SN 2004et 

NGC 6946 

5.6 ±0.3 

0.41 

0.72 ±0.03 

Over Luminous 

5,6 

SN 2005CS 

M51 

8.4 ±0.7 

0.05 

-0.10 ±0.05 

Subluminous 

7,8 

SN 2013ej 

M74 

Q q+0-4 

0.06 

1.95 ±0.06 

Normal 

This paper 


(1~) lElm hamdi et al . N200 3ll: ( 2) iLeonard et 
(5j |satiu et al.| | |2UUb[ l; (6) [A 


_ aimo^, _, 

aguire et al.| ([^ Jiu|l; (7) |P'astorelIo et al^([2i 


(3) iKrisciunas et aL 


2009|; f 4) [ Leonard e t aL]l| 2003|l 
D§lf; (8) rVmk6 el al.| < |2Ul2| | - 


time as the source cools. The bolometric luminosity de¬ 
clines very fast, by 0.4 dex in the first 30 days, and rela¬ 
tively slower by another 0.4 dex in the next 50 days. The 
UBVRIJHK luminosity is significantly different from 
the bolometric luminosity only before about -|-20 d; oth¬ 
erwise, UBVRIJHK closely resembles the bolometric 
flux. The bottom panel of Figure shows the frac¬ 
tional contribution from each of the UV, optical, and 
near-IR regions to the bolometric flux. The UV portion 
of the total flux drops from about 38% at -|-8 d to be¬ 
low 10% by -1-22 d. After this, the optical contribution 
drops very slowly and remains above 60% until the end of 
the plateau, dropping slightly during the tail phase. The 
near- IR flux contributes about 40% during the plateau 
phase, and remains almost constant in the tail. 

The log of the ratio of UBVRIJHK luminosity to 
UBVRI luminosity in the calibration sample (Table 1^ 
shows a tight correlation with B — V color. We nt 
^ogiQ{LijBVRiJHK/L ubvri) versus B — V color with a 
straight line using three of the four objects in this sample 
(Fig 14). SN 2004et is an energetic, atypical SN IIP with 
largely uncertain E{B — V)\ it is clearly an outlier, so we 
do not include it in the final fit given by Eq. below: 


logio(lT^) = (0.0856 ±0.0012) 

± (0.1056 ± 0.0012) X (R-U). (2) 

The ratios of the UBVRIJHK and calibrated 
UtoK luminosities are shown in Figure while the 
UBVRIJHK light curves are overlaid ^ calibrated 
UtoK light curves using Eq. in Figure |16| We can 
now combine the two-fold calibration: (1) llOTSE to 
UBVRI using the fit shown in Figure [T^ (fit parameters 
are given in Table [^, which gives Utol, and (2) Utol ob¬ 
tained in the first step (equivalent to UBVRI) to UtoK 
using Eq. This yields the luminosity of SNe IIP that 
have B — V and unfiltered photometry. The relative un¬ 
certainty from this procedure is about 6% or less. We 
perform the same analysis for KAIT unfiltered photom¬ 
etry. Eigure [TT] shows the final calibrated UtoK light 
curves from ROTSE and KAIT unfiltered photometry 
for SN 2013ej. Lower panels show the Ict uncertainty 
from the calibration. Although the calibration was es¬ 
tablished by limiting to data before ±102 d, it appears 
to provide reasonable estimates for the bolometric lumi¬ 
nosity even during the early nebular phase (see Fig. 

5. SPECTROSCOPY 
5.1. Key spectral features 

We present 17 optical spectra of SN 2013ej from the 
HET, Kast, and DEIMOS spectrographs in Figurej^ All 
of the spectra are corrected for the recession of the host 



B-V mag 


Fig. 14.— Linear behavior of the log of the ratio of flux with 
B-V. A tight correlation of \og^o{U BY RUHR/UBVRI) with 
B — V is seen for SNe 1999em, 2005cs, and 2013ej. The fit has 
X^/dof = 1.51. The atypical SN IIP 2004et is shown and not 
included in the fit. The derived fit is given by Eq. 


Truth versus calibration: logio(LuBVRiJHK/'-ytoK} 
20 




mean=0.0002; rms=0.0045 








-0.06 -0.05 -0.04 -0.03 -0.02 -0.01 
SN 2004et 


Fig. 15.— Histograms obtained after subtracting 
^^Bloi^UtoK/^UBVRl) from ^OgiQ{LiJBVRIJHK/^UBVRl)- 
The mean and the RMS of the difference confirm 1-2% precision 
for SN 2013ej, SN 1999em, and SN 2005cs. Applying the same 
calibration to SN 2004et (not included in the fit) yields ~ 6% 
precision. Note that the values given in the figures are in log;i^Q. 


galaxy using z = 0.002192 (NED/IPAC Extragalactic 
Databas^^ . This is consistent with that determined by 

https://ned.ipac.caltech.edu/ 
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Fig. 16.— Comparison of calibrated luminosity with measured 
luminosity. Solid lines represent UtoK light curves while empty 
and solid circles represent UBVRI and UBVRIJHK light curves. 
While 1-2% precision in the calibration is obtained for SN 2013ej, 
SN 2005cs, and SN 1999em, the outlier SN 2004et is also calibrated 
with about 6% precision. 


the Supernova Identification software (SNIDjp® Blondin 
& Tonry 20071 with a fitted redshift of O.OtS!^ Early- 


time spectra at +8 d, +9 d, and +11 d are primarily 
blue continuum, with a few P Cygni profiles of neutral 
H Balmer lines and He I lines. The opacity for all other 
ions is too low to be conspicuously observed as spectral 
features at these early phases. H I lines (Ha A6562.85, 
H/? A4861.36, and Hy A4340.49) are very broad at early 
times. The strength of emission component of H I lines 
decreases with time. At +11 d, O I A7775 is glimpsed. A 
week later at +19 d, several strong absorption signatures 
of SNe H appear. 

Interestingly, the absorption line to the blue of Ha is 
unusually strong. This feature, which we identify as Si H 
A6355 (see Section 5.3), subsequently becomes stronger 
until +19 d in our sample. It appears as a small absorp¬ 
tion notch a t +44 d and disappears by +48 d. |Valenti| 
et al. (2014) showed this feature to get stronger tnan Ha 
until -j- 21 d in their dataset, and to become weaker than 
H a by +23 d. The Si H identification was also favored 


et aT 
in other 


by Valenti et a l. (2014),|Bose et al. (2015a), and Huang 
(|2015l). Si 11 has seemed to occur much later 


SiNe IIP. While this strong early appearance of 
Si H has not been observed previously, it may have been 
marginally detected at +10 d and +12 d and was not ob - 


served after +25 d in SN 2006bp (Quimby et al. 


2007). 


Si H is comparatively much stronger than in Si\ zOOhbp 
at similar epochs. For SN 2006bp, the Si H velocity pro¬ 
file evolves faster than that of Ha before +25 d, whereas 
for SN 2013ej, it is more smooth and evolves more slowly 
than the Ha velocity. All these factors make SN 2013ej 
exhibit unusual and strong early Si H. 

As the ejecta expand, the subsequent spectral evolu- 


http://people.lam.fr/blondin.stephane/software/snid/ 


tion of SN 2013ej shows typical SN HP singly ionized 
lines of Ca H, Fe H, Ti H, Sc H, and Ba H. The He I 
A5876 line gradually gets weaker until +15 d and is not 
seen at +19 d. This is evidence of the temperature 
decreasing below the critical temperature of excitation. 
More iron-group elements start to appear, corresponding 
to the commencement of the plateau phase where the 
photosphere penetrates deeper into the envelope. The 
same disappearance of He I at around + 16 d was also 


seen in SN 1999em (Leonard et al. 2002 

) and recently 

in SN 2012aw ( 

Bose et al. 

20T3 

). 'i'he Na 1 U lines 


are considerably stronger than the lines of other neu- 


tral elements, 

presumably coming from non-LTE effects 

(Hatano et al. 

1999 

). The Na I D feature is observed 

after +19 d anc 

1 is probably blended with He i at +15 d. 


We do not observe any narrow lines of Na I. No obvi¬ 
ous evidence of high-velocity features (HVFs) is seen in 
our spectra. These observations may indicate negligible 
interaction of the ejecta with the circumstellar material 
(CSM). After +19 d, the Ca H near-IR triplet can be dis¬ 
sociated to at least a doublet at 8520 A and a singlet at 
8662 A; however, the profile is well blended before +15 d, 
and we adopt this as a single Ca H near-IR profile to de¬ 
termine the change in ion velocity with time. 


5.2. Spectral Homogeneity in the UV 

SNe HP are known to exhibit a remarkable h omogene¬ 
ity in their U V spectra, as first pointed out by Gal-Yam 
et al. (2008). They found that the early-phase UV spec- 
tra (2000-3000 A) of SNe 1999em, 2005ay, and 2005cs 
are very similar, both in the shape of the continuum 
as w ell as in the visi b le spe ctral features. In compar¬ 


ison, 


Ben-Ami et al. 


SNe Hb, which are t 


(2015) recently pointed out that 
lought to have thinner H-rich en¬ 


velopes than regular SNe HP, display relatively strong 
diversity in their UV spectra. 

The paucity of well-observed SNe HP having early-time 
UV spectra impedes an in-depth study of this homogene¬ 
ity verses diversity issue at present. It is therefore im¬ 
portant to increase the size of the early-time UV sample. 
SN 2013ej is a valuable addition to this sample because 
of its relative proximity, which enabled Swift to obtain 
near-UV spectra with its UVOT/UGRISM instrument 
(see Fig.[^. Figure flSl compares the +8 d and the +11 d 
spectra to those of other SNe H taken at similar phases. 
All of these spectra are corrected for interstellar extinc¬ 
tion and scaled to match the fluxes in the region 2500- 
3000 A. 


Figure 18 reveals that SN 2013ej nicely fits into the 
frameworkof the UV spectral homogeneity of SNe HP, 
at least around 10-12 days after explosion. We find that 
the similarity is not evident for spectra taken at ~ 1 
week after explosion (Fig. 18 top panel) in our sample. 
Both SN 1987A and SN 2(jD5cs showed some differences 
with respect to the spectrum of SN 2013ej at this phase, 
although the rise of the UV flux in the SN 2005cs spec¬ 
trum below 2500 A may not be real. Close inspection of 
the UVOT/UGRISM frames revealed that this spectrum 
was contaminated by emission from the zeroth order of 
a nearby source. Moreover, SN 1987A, which shows a 
sharp cutoff in the UV flux below 3000 A, was not a 
typical SN HP, as it had a blue supergiant progenitor. 
Nevertheless, the spectra taken around 11 + 1 days after 
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Fig. 17.— UtoK luminosity obtained from calibrating broadband and unfiltered photometry from KAIT, Konkoly, and ROTSE data for 
SN 2013ej. Data points UBVRIJHK in the top panel are integrated U from Swift, KAIT BVRI, RATIR JH, and an estimated K flux 
(see text). Shaded regions in the residual plots indicate Icr uncertainty from the RMS of the residuals. For KAIT data, we have added 
systematic uncertainty based on the differences of KAIT UBVRIJHK and UtoK values derived from UBVRI luminosity using Eq. 
which is derived using Konkoly data. The highest offset at very early times is a result of the high UV contribution to the total flux. 


- 

ROTSE Unf. Calibrated UfoK"-Konkoly UBVRIJHK - 
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explosion confirms the o bse rved similarity nicely (Fig. |18[ 
bottom panel). Figure 18 also illustrates a SN Ilb 


spectrum at a similar epoch; it differs significantly from 
the SN IIP sample. 

5.3. Line Identification and Spectrum Modelling 
Line identifications of most of the features in Section 


5.1 were first driven by the study of Hatano et al. (19991 
on ion signatures of SN spectra. Addit ional study and 
confir mation was performed by Syn++ (Thomas et al. 
20111 modelling of a few of the optical spectra as shown 
in Figure While the synthetic modelling produces 
most of the ionic signatures, the most obvious Ha profile 
is not reproduced. This reflects the limitation of a purely 
scattering code: the emission is underestimated because 
it does not account for the emission due to recombina¬ 
tion cascades. For Ha, there may also be a significant ef¬ 
fect from non-thermodynamical equilibrium (NLTE) and 
time varying effects. The absorption notch blueward of 
the Ha line in the -1-11 d and -1-19 d spectra is fitted with 
the Si H line, and we have obtained the fit as shown in 


Figure]^ One could argue that this feature is an HVF 
of Ha, out then we would expect to also see HVFs of 
other Balmer lines. Taking an HVF input with such a 
high velocity, we were unable to reproduce a decent over¬ 
all fit. Because of the lack of an HVF for other Balmer 
lines and no HVFs seen in t he near-IR spectra, a s re¬ 
ported at similar epochs by Valenti et al. (2014), the 
HVF hypothesis is disfavored. Clearly, the Sill identifi- 
cation hypothesis will be settled with high er confidence 
only fro m more realistic modelling. While |Bose et al. | 
(2015a) also identified the blueward notch as a Si 11 fea¬ 
ture, tney have incorporated HVFs for H I lines, albeit 
blended with the photospheric component, accounting 
for the broad Balmer lines beyond -1-42 d in their sam¬ 
ple. After -1-15 d, lines of intermediate-mass elements 
and iron-group elements start to appear. The s-process 
products Ba H and Sc H are also seen from -1-19 d until 
the last spectra in the plateau at -1-94 d. 


5.4. Velocity Evolution 
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Rest wavelength (A) 



Rest wavelength (A) 


Fig. 18.— Evidence of spectral homogeneity of SNe IIP in the UV 
at early times. Top panel: SN 2013ej UV spectrum compared with 
the atypical SN 1987A, which shows a sharp UV cutoff, while the 
excess flux of SN 2005cs below 2500 A is suspected to be coming 
from a different source. Bottom panel: Homogeneous SN IIP UV 
sample at ~ 12 d. Also shown for comparison is a SN Ilb spectrum, 
which is clearly distinct from the rest. 


While the average ejecta velocity is a direct tracer of 
kinematic properties, the photospheric velocity not only 
provides compositional clues but also traces the size of 
the photosphere, thereby aiding distance measurements 
(e.g., EPM). The photospheric velocities at early spectral 
epochs are estimated by fitting the He I A5876 feature. 
After +19 d, the Fe 11 A5169 line is most indicative of the 
photospheric velocity, since the minimum of the absorp¬ 
tion profile tends to form near the photosphere (IBranch 
et al. |[200^ . 


The velocity evolution of some of the strongest ions is 
presented in Figure Each line absorption feature is 
fitted by a Gaussian profile, and the minimum is con¬ 
verted to velocity using the relativistic Doppler equa¬ 
tion. The Ha line is decelerating more slowly than H/3 
and other metallic ions as expected, but the H I lines 
sh ow a flat velocity pr o file, which was also poi nted out 
by Bose et al. (2015a). Poznanski et al. I (20101 demon- 
strated a correlation of velocity of the Fe 11 A5169 line 
{vpe II ) with that of the H/3 line (uh^a) using 28 optical 
spectra of 13 SNe HP c overing 5-40 days after explosion. 


Takats & Vinko (2012) have extended the validity of this 
relation to phases beyond 40 days. We have examined 
this behavior by taking four optical spectra of SN 2013ej 
from 15-48 days after explosion. Analysis of Fe 11 is not 
justified before +15 d in our sample. We find that the 
velocities from SN 2013ej spectra are consistent with this 


TABLE 7 

Photospheric velocities of SN 2013ej determined by 
Syn++ fitting. Phases are rounded to the nearest day 
SINCE explosion 


MJD 

Phase 

’^phot 

Uncertainty 


(days) 

(km 8“-'-) 

(km s“^) 

56505.5 

+8 

10200 

1000 

56506.5 

+9 

9700 

1000 

56508.5 

+11 

8800 

800 

56516.5 

+19 

7660 

600 

56541.5 

+44 

4700 

500 

56545.5 

+48 

4900 

400 

56566.5 

+69 

3200 

400 

56570.5 

+73 

3740 

500 


correlation, as can be seen in the right panel of Figure!^ 
We obtain a linear relation of upe ii = 0.85+0.03 uh/I iot 
SN 2013ej, in agreement with vpp n = 0.84 + 0.05 uh /3 
as obtained by Poznanski et al. (2010). 

In order to rehne the photospheric velocity using more 
than one feature, we first approximate the velocity by es¬ 
timating from the absorption minimum of the He 1 A5876 
line for the two earliest spectra and the Fe 11 A5169 line 
for the later spectra. We then performed synthetic spec¬ 
tral modelling with Syn-|—I- and extracted photospheric 
velocities from the model. Velocities obtained from the 
fits are given in Table 


6. DISTANCE DETERMINATION 
Recently, the distance to the SN 2013ej host, M74, was 


of D « 7 + 2 Mpc (Sharina et al.|1996 

Vinko et al. 12004 

Van Dyk et al. 

20U6D to D « 9.5 + ( 

1.5 Mpc ( 

Zasov & 

Bizyaev||199t) U 

livares E. et al.||2010) 

; see Tab 

e|^for a 


have been applied for various samples of SNe HP (Hamuy 

et al. 

2001 

Leonard et al.|20021|Dessart et al.|2U081|Jones 

et al. 

2009 

Vinko et al.|2012||Bose & Kumar|2014|[Takats 

et al. 

2014 

). We present the application of the version 


from two SNe that occurred in the same host galaxy, 
claiming the uncertainties of EPM can be reduced and 
the reliability of the derived distance improved. Thus, 
we take the advantage of having the necess ary data for 


both SN 2013ej (this paper) and SN 2002ap (Vinko et al. 


2004), although the latter object is a broaddined SN Ic 


for which the application of EPM may not be fully jus¬ 
tified. Despite the complications arising in modeling the 
atmospheres of such stripped-envelope (SE) CC SNe, we 
show below that the combination of the two datasets re¬ 
sults in surprisingly consistent results, and the inferred 
distance is in very good agreement with recently pub¬ 
lished independent estimates. _ 

Foll owing the procedure described by [Vinko et al.| 


(2012), the basic equation for EPM is 


t = D X 


(—) 

\ '^phot / 


to, 


( 3 ) 


where t is the time, D is the distance, 0 is the angular 
radius of the photosphere, Uphot is the velocity of the 
photosphere at t, and to is the moment of shock breakout. 
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Fig. 19.— Example SynH—h fits of SN 2013ej spectra are shown in blue while data are in black. The fits mostly reproduce the observed 
features. The inability to accurately reproduce the H I line profile is perhaps a limitation of the model being purely scattering based and 
not accounting for the emission from recombination cascades and the NLTE effects. SN 2013ej exhibits most previously identified SN IIP 
spectral features. 
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Fig. 20.— Left Panel : SN 2013ej velocity evolution of strong ions. Empty black ci rcles are the photosp heric velocity derived from 
Syn+H- fits. The dashed line is a fit to the photospheric velo city using the method of IVinko et al.l (|2012|| (see text). Right Panel : 
Demonstration of correlation of upe n with as suggested by|Poznanski et al.| (|2010||. I'he stiaded regions indicate the Ifj region of the 
correlation: gray, [Poznanski et al.| l [2010| |; navy, SN 2013ej (this paper). 
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We estimate 9 from the bolometric light curve by using 


1 


CiT) 


I /bol 


(4) 


where C(^) is the dilution factor describing the alteration 
of the pure blackbody flux in a scattering- dominated 
SN atmosphere as a function of temperature (Eastman, 


Schmidt & Kirshner 1 1996 Dessart & Hillier 2UU5|) and 

/bol is the apparent bolonietric flux. I 'or SN~^2U13ej, we 

u sed th e dilution factors determined by [Dessart fc HilliCT 


[ (2005|), which are valid for H-rich SNe IIP, but not for 


the M-free SE SN 2002ap. Since the atmospheres of such 
SE SNe are much less known, we set C = 1 as a first ap¬ 
proximation. Note that the usage of C = 1 worked sur¬ 
prisingly well when calculating th e distance to the Type 
Ilb SN 201 Idh (Vinkoet al. 2012). Since the ejecta of the 
Type Ic SN 2002ap contained practically no H, unlike the 
Type lib SN 201 Idh, the dilution of the blackbody flux 
due to Thompson scattering on free electrons might be 
even less strong than in the case of Type II SNe 2013ej 
or 201 Idh. Thus, setting « 1 may be a physically 
realistic approximation for SN 2002ap, although its full 
justification would involve the computation of an NLTE 
model atmosphere for SN 2002ap which is beyond the 
scope of this paper. Nevertheless, we estimate the prob¬ 
able amount of the systematic error of the distance in¬ 
troduced by the assumption of ^ = 1 below. 

The estimates of 9 were based on the bolo metr ic light 
curve of SN 2013ej, as described in Section [431 More¬ 
over, we applied the /boi fluxes of SN 2002ap similarly, 
afte r coinbining the opt i cal lig ht c urves from [Foley et| 
aTl (|2003|) , jPandey et al.j ([2003 1, an dIVinko et al. | (|2004 |) 
with the near-lK measurements by Yoshii et al. (2003). 
In the latter case, the UV contribution was estimated by 
assuming zero flux at 3000 A and a simple linear SED 
between 3000 A and the U band. This approximation 
was justified by the shape of the spectra of SN 200 2ap as 
they declin ed below 4500 A toward the blue (e.g., 
er^[2004| . 


Vinko 


i'he application of Eqj^ and [5 requires Uphot and 
values at several epochs, typicmly during the first SO¬ 
SO days after explosion. These can also be estimated 
directly from the observations. In case of SNe Ic like 
SN 2002ap, the applicability of EPM is limited to no 
longer than a few weeks. 

F or SN 2013ej, values of Tes were obtained in Section 
|4.4| For SN 2002ap, we applied the re ddening value of 
E{l3—V) = 0.09 mag ( Vinkd et al.|2004 1 and the relation 

(5) 


Teff = - 0.122(5-P) -h 3.875 


based on the SN Ic-BL models by Mazzali et al. (2000). 

In order to increase the sampling ot the velocity curve 
of SN 2013ej, we fit the velo city curve for SNe IIP derived 
by TaMts & Vinko (2012) to the velocities obtained in 
Table [Tj 'i'his method involves velocity modelling as a 
power-daw expansion of phase, given a model velocity at 
some epoch. This model velocity is generally derived 
from synthetic modelling of the observed spectra or by 
direct measurement from the absorption pro file of lines 
like Fe II A5169. See Takats & Vinko (2012) for a more 
detailed discussion of this kind ot velocity measurement 
in SN IIP atmospheres. The result of this fitting was 


TABLE 8 

Physical Quantities Derived for EPM. 


time 

e 

^/■^phot 

Uncertainty 

(days) 

(10® km Mpc“^) 

(day Mpc“^) 

(day Mpc“^) 

SJN 2013ej 

8.60 

9.44 

1.08 

0.12 

10.60 

10.10 

1.24 

0.12 

13.60 

11.45 

1.52 

0.17 

14.60 

12.78 

1.74 

0.19 

15.60 

12.98 

1.81 

0.20 

16.60 

13.30 

1.90 

0.20 

19.60 

13.93 

2.13 

0.22 

20.60 

15.31 

2.39 

0.24 

24.60 

16.39 

2.77 

0.28 

25.60 

15.80 

2.72 

0.28 

SN 2002ap 

4.89 

11.96 

0.44 

0.34 

6.48 

12.34 

0.67 

0.30 

7.48 

13.00 

0.78 

0.32 

9.87 

14.70 

1.06 

0.35 

10.87 

15.17 

1.17 

0.37 

11.27 

15.53 

1.19 

0.38 

12.87 

16.15 

1.55 

0.40 

13.47 

16.11 

1.51 

0.39 

13.86 

16.10 

1.69 

0.39 


applied in the procedure of EPM. 

For SN 2002ap, we adopted t he velocities based on the 
Si II A6355 feature as given by [Vinko et al.j ([2004[) . 

The derived quantities needed tor EPM are shown in 
Table [^ The moments of the explosion were set to be 
to = M JD 56496.9 (2013 July 23.9 UT), as derived in 
and to = MJD 52302.0 (2002 Jan 28.0 UT) 
3ej and 2002ap, respectively. 


Section 4.1 
for SNe^ 

The ht of Eq. [^ to the data in Table [^ was performed 
assuming 9/v as the independent variable, with either 
keeping to fixed at the values given above or letting it 
float. The first fit resulted in 5 = 8.86 ± 0.21 Mpc, 
while the second one gave D = 9.09 ± 0.30 Mpc with 
Ato = —0.59 ± 0.47 being consistent with our estimated 
to. Alternatively, choosing t as the independent variable, 
one may get D = 8.93 ± 0.10 Mpc and D — 9.25 ± 0.30 
Mpc with Ato = 0-09 ± 0.48 days for fixed and float¬ 
ing to, respectively. The weighted average of these four 
values gives D = 8.96 ± 0.08 Mpc. Accounting for any 
systematic effect that could have been introduced from 
our derived to in Section |4.1[ we take lower and upper 
bounds for to as —1.3 d (obtained from floating index, 
see Section 4.1) and -1-0.9 d (Lulin detection epoch) from 
the derived to. Fitting with adjusted lower and upper 
bounds of to, we get 0.35 Mpc higher distance and 0.60 
Mpc lower distance, respectively. We adopt these offsets 
as the systematic uncertainty and add in quadrature with 
the statistical uncertainty obtained above. Thus, O.Oj'jQ g 
Mpc is adopted as the final distance estimate of M74 
from EPM using two SNe. This derived distance can be 
found from the inverse of the slope of the line shown in 
Figure [^ 

In order to test the effect of choosing C = 1 artificially 
for SN 2002ap, we repeated the fitting process described 


above after ap plying the dilution factors of [Dessart fcj 
Hillier ( 2005( ) to the SN 2002ap data as weE This is 


clearly an overestimate of the effect of electron scatter¬ 
ing (i.e. an underestimate of the C values) in a Type Ic 
SN atmosphere, which may be somewhat less scattering- 
dominated than a H-rich Type HP atmosphere. How¬ 
ever the amount of the systematic error introduced by 
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TABLE 9 

Recent Distance Estimates eor M74. 


Fig. 21.— Distance measurement of M74 using SN 2013ej and 
SN 2002ap. The black solid line shows the final result that yields 
a distance of 9.0 Mpc, while the dotted lines mark its uncertainty. 


such a strong dilution might be useful for constraining 
the real uncertainty of the distance due to the approx¬ 
imate dilution factors. Assigning the C{T) values from 


(plotted with red triang. 


Dessart & Hillier (2005) to SN 2002ap move those data 


es in Fig. [^by about ~ Icr up¬ 
ward, reducing the consistency between the two datasets. 
Performing the same fitting process as described above, 
we get D = 8.3 ± 0.6 Mpc, i.e. less than 2a differ¬ 
ence from the previous distance estimate from ^ = 1. 
Since this test uses potentially underestimated values of 
( for SN 2002ap, we conclude that the uncertainty caused 
by the inaccurate knowledge of the dilution factors for 
SN 2002ap probably does not exceed the ±0.5 Mpc un¬ 
certainty estimated above. _ 

Using an in dependent dataset for SN 2013ej only, |Rici>| 
mond (2014) applied the standard version of EPM to 
get U = 9.1 ±0.4 Mpc. From the values given in their 
Table 5, the uncertainty of that distance appears to be 
around ~ 0.8 Mpc instead of ±0.4 Mpc as noted. In ei¬ 
ther case, this is in very good agreement with our result. 
Furthermore, we estimated the distance using the bolo- 
metric calibration for ROTSE unfiltered fluxes derived in 
Section [±^ For this case, we include only data points be¬ 
yond ±15 d for SN 2013ej, as earlier data would include 
significant flux below the U band that is not included in 
the calibration procedure. We get an EPM distance of 
9.7 ± 0.6 Mpc, using the calibrated fluxes from ROTSE 
for SN 2013ej, in agreement with the previous deriva¬ 
tion. Note that we have not combined SN 2002ap values 
in this case. Allowing upper and lower bounds to to as 
before, we get the final distance estimate from calibration 
of ROTSE data to be 9.7^° ® Mpc, which is consistent 
with our preferred 9.0to(6 result. 

7. EXPLOSION PROPERTIES 
7.1. Ni Mass 

The end of the plateau phase is believed to indicate 
the full recombination of hydrogen when the ionization 
front, and thus the photosphere, reaches the bottom of 
the hydrogen envelope in the ejecta. After this epoch, 
the light curve proceeds into a nebular phase. The sub¬ 
sequent luminosity is driven by the radioactive decay of 
elements that were produced during the explosion, so 
the light curve shows a characteristic exponential decay 


Method 

D 

Heterence 


(Mpc) 



T-F 

9.68 ±1.63 

ri-ully 

J11988II 


BBSG 

7.31 ± 1.23 

Ibnanna e 

t ±1.111996 


Disk gravitational 

9.4 

Ziasov 6c tsizyaev UiyiJ 

Ibl 


stability 
Light echo 
SCM 

(SN 2003gd) 


7.2 

9.91 ± 1.2 


IVan Dvk et al. > 

2006b 

|Ulivares iii. et a . 

l|2U10|) 


EPM (SN 2002ap) 

6.7 

JVinko et al. 

0 

0 

TRGB 

10.2 ± 0.6 

fjang & Uee 

12U14| 


EPM 
EPM: ROTSE 
calibrated 


9.03 

9 7 +°'^ 
' -0.7 


present paper 
present paper 


of flux output. This suggests that the gamma-rays and 
positrons from radioactive decay of ®®Co thermalize in 
the ejecta. Here, we first assume the full trapping of 
gamma-rays and positrons in the ejecta. As the mass of 
freshly synthesized Ni should be proportional to the tail 
luminosity, we use the findings from the literature and 
use scaling to determine the nickel mass (Mni) for SN 
2 013ei. 


Bose et al. (2013) derived the Mn; for SN 2012aw 


using the UBVRl fight-curve tail luminosity. As we 
will see below that the decline rate changes after ±183d, 
we linearly fit the UBVRl light curve from ±120 d to 
±183 d and extrapolate to find the luminosity at 240 d to 
make a direct comparison with their result for SN2012aw. 
The luminosity, L(240 d) for SN 2013ej, is estimated to 
be 1.32 ± 0.05 x 10"^° erg s“^, while that for SN 2012aw 
was found to be 4.53 ± 0.11 x 10^° erg s“^. The ratio is 
calculated to be 0.29 ± 0 .02. Noting Mni fo r SN 2012aw 
to be 0.058 ± 0.002 Mq (|Bose et al. |2013|), we calculate 
for SN 2013ej, Mni = 0.017 ± 0.001 Mf^. Alternatively, 
we use the method of Hamuy (2003) to calculate the Ni 
mass from the tail luminosity \Lt), using the equation 


Mni = 


7.866x10 "‘"^Ltexp 


[it - fo)/(l + z) - 6.1 


111.26 


Mq. 

( 6 ) 

We calculated Lt at 20 epochs from ±120 d to ±183 d 
using the late-time U-band magnitude. We ad opt the 
s ame b olometric correction of 0.26 mag from |Hamuy| 
[] (2003). The weighted mean tail luminosity is calcu- 
lated to be 5.82 ± 0.26 x 10'^° erg s“^, corresponding 
to ±157 d. The Mni value is then calculated to be 
0.018 ± 0.002 Mq. Following the same procedure, but 
using the bolometric light curve derived in Section |4.5[ 
we get Mni to be 0.019 ± 0.003 M©. We take the 
weighted mean of the above three results as our best 
estimate of the synthesized radioactive material. This 
y ields Mni =0.0 18 ± 0.001 Mq 


Hamuy (2003) used a large sample of SNe HP to study 
the correlation between the Ni mass and mid-plateau 
(±50 d) photospher ic v elocity. From the modelling as 
described in Section [5^ we have derived a photospheric 
velocity of 4500 km s“^ at ±50 d. Our results for A/n , 
and V 50 are consistent with the results of Hamuy (|2003 ). 

We note, however, that from the late time data of 
SN 2013ej from the KAIT and Nickel telescopes, we ob¬ 


serve two distinct slopes in the tail (see Fig. 22). To 


estimate the time of slope break, we fit the late time 
bolometric flux beyond ±120 d with a broken exponen- 
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Fig. 22.— Fit of the bolometric light curve of SN 2013ej beyond 
+120 d with a broken exponential law. The solid line is the fitted 
model while the dashed lines are the extrapolation of each expo¬ 
nential models in either direction. The crossover is at break time, 
+183.4+ 15.7 d. 


tial law of the form, 

F{t) = SAe~^i 


l + e“(‘-*'”-) ^ (7) 


where, ti and T 2 are characterstic times for the Hrst and 
second exponential prohles, A is the initial flux, is the 
break time, a is the smoothing parameter and S is the 
scaling factor, given by. 


S = (l + e-"**”-) 


t V 1-1 T2 ^ 


( 8 ) 


Eq. [^has been analogously applied to study the radial 
profi le of surface brightness fr om the disks of galaxies 
(e.g.|Munoz-Mateos et al.|2013). The best fit parameters 
are f+nd to be, t, = 73.89±5.00 days, T 2 = 94.73± 1.39 
days, tbr = 183.28 ± 15.67 and a = 0.23 ± 1.14 with 
the /dof from the fit to be 1.55. The slopes before 
and after the break point are obtained to be, Slopei = 
0.015 ± 0.001 mag day^ and Slope 2 = 0.011 ± 0.001 
mag day“^. While Slopei is found to be much steeper, 
Slope 2 is closer to the ^®Co —>■ ®®Fe decay rate. 

To see the effect of these two distinct decline behaviors 
on the initial Ni mass, we further fit the late-time bolo¬ 
metric lig ht curve of SN 2013e i wit h the simple model de - 
scribed in Vinko et al. (2004) and Valenti et al. (2008). 
This model assumes an optically thin ejecta heated by 
the radioactive decay of ^®Ni and ®®Co. The decay en¬ 
ergy is emitted in the form of gamma-rays and positrons, 
which may be partially trapped in the ejecta, thermal- 
ize and emerge again as low-energy (mostly optical or 
near-infrared) photons. The deposition function for the 
gamma-rays at a given epoch t can be expressed as 

= l-exp[-{'^n (9) 


= 1-e"^^ 


where is the optical depth for gamma-rays in the whole 
ejecta. The timescale of the gamma-ray optical depth 
decrease (Wheeler et al.|[20T5l is 


To (7) = 


( 10 ) 


where Kj is the gamma-ray opacity, Mgj is the ejecta 
mass and C is a constant depending on the density dis¬ 
tribution in the ejecta. For simplicity, we assumed a 


Fig. 23.— Fit of the radioactive decay model of Eq. [D to the 
bolometric light curve tail of SN 2013ej. Pre-break time and post 
break time lightcurves are fitted separately. The upper panel shows 
the fit taking data between +120 d and the break point at +183 d. 
The lower panel shows the same fit taking data beyond the break 
time. The leakage of both the gamma-rays and positrons were 
taken into account in the model (see text). The red line represents 
full trapping of gamma-rays and positrons on both panels. 


constant density ejecta, which implies C = 9/407r. The 
deposition function for positrons, takes the same 
form, except for the opacity. For the gamma-ray opac¬ 
ity, we adopted = 0.027 cm^ g~^ and for positrons, 
we set K+ = 7 cm^ g“^ (e.g. Colgate et al.||l980 Valentif 
eral][200^ . 

With these definitions, the late-time bolometric lumi¬ 
nosity can be expressed as 


Lboi = ArNi[(S'Ni(t) -f O.925'co(0)+ 
-f {0m + 0m*D^)Sco{t)D+], 


( 11 ) 


where Mni is the initial mass of the radioactive ®®Ni 
synthesized during the explosion, 5'Ni and Sco are the 
functions for the total en ergy input from the Ni- and 


Co-decay, re s pectiv ely (see Szalai et al. (2016); Branch 
& Wheeler (2016) for further discussion^ TTis equa¬ 


tion corr ects for the typograp hical error in the expression 
given by |Valenti et al.| (|2008|) , and accounts for the par¬ 
tial trapping of both gamma-rays and positrons via the 
deposition functions given above. Since this light curve 
model assumes instantaneous release of the thermalized 
deposited energy from radioactive decay, witho ut con¬ 
siderin g any photon diffusion unlike the model of lArnetT] 
(1980), it is applicable only when the ejecta is almost 
fully transparent in the optical, i.e. during the nebular 
phase. 


The fit of Eq. 11 to the bolometric lightcurve is plot¬ 
ted in Fig. ^ found that before -1-183 d, SN 2013ej 
exhibited a steeper decline in the bolometric ligh t curve 


than th e rate of the Co decay, as also f ound by Huang 


et al. (12015 ) and Bose et al. (2015a); however, our 


extended photometry, revealing a shallower decay rate 
has an implication for both the gamma-ray opacities 
and the estimated Ni-mass. From the fit beyond the 
break point of -1-183 d, the values inferred for the ini¬ 
tial nickel mass and the gamma-opacity timescale are 
Mm = 0.013±0.001 Mq and To( 7 ) = 465±18 days. This 
timescale is s ignificantly longer t han the value of ~ 173 
day found by Bose et al. (2015a) from the data between 
-1-100 and -|-2(j0 days. (Jur longer To( 7 ) suggests that the 




















































20 


Dhungana et al. 


light curve of SN 2013ej was probably not yet settled on 
the radioactive tail before +183 days. The Ni mass de¬ 
rived from the Ht before break time, Mni = 0.019 + 0.001 
Mq is consistent with our calc ulations above u sing meth¬ 
ods of|Bose et al. | (|2013|) and |Hamuy | (|2003D, and with 
previous independent estimates by Huang et al. (2015) 
and Bose et al.(2015), both assuming full gamma-ray and 
positron trapping in about the same time window, but all 
those estimates are probably overestimates and the lower 
value derived by restricting the analysis beyond +183 d 
is more likely to be correct. 

7.2. Explosion Physical Parameters 

To determine accurate explosion properties, a detailed 
hydrodynam ic study is required. Recently, [Huang et 


aj^ Jl2 015 ) have done such a study while Bose^^aL 


use a semi-analytic approach following | Arnett^^ 


h'u 11989). In order to make an approximate estimate , 
we use the approach of Litvinova & Nadyozhin (1985). 


Even though there are complications concerning radioac- 
tive heating effects, inclusion of higher explosion energies 
in the models, simple physical assumptions will still be 
useful to compare the explosion parameters with those 
of the more elaborate studies. 

From plateau length (At), the mid-plateau absolute V 
ma gnitude My, and the correspond ing expansion veloc¬ 
ity, Litvinova & Nadyozhin ( |1985[ ) derive the explosion 


energy, ejected mass, and pre-SJN radius. For the plateau 
midpoint of SN 2013ej, we use (tpeak + ip)/2; where tpeak 
is the epoch of peak brightness in V, estimated to be 
+15 d using Gaussian Process regression of V band data 
until +30d, and tp is the epoch at which the plateau 
ends. It is nontrivial to precisely locate the end of the 
plateau. To determine this epoch, we again perform a 
Gaussian Process regression on the bolometric lightcurve 
from +80 d to +130 d, and determine the point of inflec¬ 
tion from the obtained ht to be +109 d. We consider this 
to be the end of the plateau, tp. We therefore estimate a 
plateau of length 94+7 with -1-62 d as the midpoint. The 
value of Uph at +62 d is found to be 3800 + 500 km s“^ 
from the ht described in Section [5^ while My is deter¬ 
mined to be —16.47 + 0.04 mag from linear interpolation 
of V magnitudes from +50 d to +70 d. Using these val- 


ergs, and a pre-SN radius of 250 + 70 Rr7^ 

based 

on Litvi- 

nova & Nadyozhin 

1_^ ___ 

(1985). 

Bose et al. 

(2015a 

C 

1 hnd an 


of 450 R0. Huang et al. (2015) obtain values ranging 


over 0.7 — 2.1 x 10°^ ergs and 230 - 60 0 Rp, respec¬ 
tively. Given the unstated uncertainty in |Bose et ah 


(2015aI, and the range from Huang et al. (|2015|), our 


calculations appear to be consistent with theirs. Uur ex¬ 
plosion energy and pre-SN radius are also in agreement 


Hamuy 

(2003 

) and Nadyozhin (2003 

, both of which 

use the same 

jitvinova &: Nadyozhin 

(1985 

) relations 


We further use these relations to determine the ejecta 


mass, M fi, for SN 2013ei to b e 13.8 + 4.2 Mq. Hamuy 
[ (2003) and Nadyozhin (2003) obtain Mgj in the range 


14-56 Mq and 10-30 Mq, repsectively, which encom¬ 
passes our result. As for the explosion energy and pre- 
SN radius, the ejecta mass is highly sensitive to plateau 


length; better knowledge of the plateau length c an yield 
more accurat e res ults. For SN 2013ei specihcally, |Huang| 
et al. (2015) and Bose et al. (2015a) hn d ejecta masses 
of 10.6 Mq and 12 + 3 Mq, respectively. |Nagy & Vink6| 
(20161 have derived an ejecta mass of 1(1.6 Mq from light 
curve modeling using a two-component model incorpo¬ 
rating a dense inner core and an extended low mass en¬ 
velope. These three esti mates are all consistent w i th our 
measurement based on Litvinova & Nadyozhin (19851 
relations. If we assume a remnant mass of 1.4 Mq, our 
measurement for the hnal pre-explosion progenitor mass 
is 15.2 + 4.2 Mq. 

Other measurements have been performed of SNe HP 
progenitor mass in general or SN 2013ej specihcally. Gare 
is required in comparing progenitor masses, however, 
as some correspond to initial pre-explosion masses or 
ZAMS masses rather than the hn al pre-explosion pro - 
genitor mass we calculated above. Fraser et al. (2014), 
for instance, found the ZAMS mass to be in the t- 
15.5 Mq range for SN 2013ej, consistent with an M- 
type supergiant, fr om the archival HST imag es. Using 
X-ray observations, Chakraborti et al. (2015) derived a 
ZAMS mass of 14 Mq , accounting for the derived steady 
mass loss of 3 x 10“° Mq yr~^ over the last 400 years. 
Given the uncertainties, these measurements are consis¬ 
tent with our hnal progenitor mass of 15.2 + 4.2 Mq. 


For th e general population of SNe HP, Smart! et al. 
(2009) obtain a ZAMS mass range of 8-17 Mq. Use 


of nebular phase modelling (e.g. jJerkstrand et al. [2014 1 
can also independently provide tighter constraints on the 
ZAMS masses of these events. However studies involv- 


ing hydrod ynamical modelling (e.g. Utrobin & Chugai 
2009[ 2013 ) and stellar evolutionary models (e.g. Smart! 
et al. ^2l)09[ ), along with nebular spectra modellirig, have 
shown conhicts in the derived initial mass of the progen¬ 
itor stars. 

8 . DISCUSSION AND CONCLUSIONS 

We present extensive photometry of the nearby 
SN 2013ej at UV, optical, and near-IR wavelengths. We 
also discuss well-sampled UV and optical spectroscopy 
from +8 d to +135 d after explosion. SN 2013ej looks 
kinematically similar to other normal SNe HP, but it 
also exhibits some unique features compared to a broader 
sample of SN HP, such as a steep plateau, early ap¬ 
pearance of strong Si 11 A6355, and a hat Ha velocity 
prohle. Such features hint at an intermediate class be¬ 
tween SNe HL and HP, or probably a continuum in the 
distrib ution of these CC SNe. From a large sample of 
SNe H, Anderson et al. (20141 did not hnd any evidence 
of bimodality in the distributions of many photometric 
properties they studied based on the V band, thereby 
suggesting a continuum in the properties of SNe H. 

SNe HP show a wide range of plateau duration. Even 
the typical SNe H P that exhibit hat p lateaus, like 


SN 2006 bp (<73 d, Quimby et al. 2007), SN 2013ab 
(^ 80 d, iBose et al . 2015b|), and SN 2003hn (^ 75 d, 
Bersten et al. |20ll|) , nave shorter plateau duration com- 


pared to that of SN 2013ej. This indicates that the en¬ 
velope mass of SN 2013ej is not atypical. The amount 
of decrease of luminosity from the end of the plateau to 
the radioactive phase may be related to production of 
Ni in the ejecta, but the value of Mni cal culated for SN 
2013ej is also not atypical of SNe HP (e.g., Hamuy [2003 
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TABLE 10 

Calculated Physical Parameters of SN 2013ej. 


Parameter 

Bose et al. {2015a) 

Huang et al. (2015) 

Fraser et al. (2014)“ 

Valenti et al. (2014) This paper 

Explosion Energy (10^^ ergs) 

2.3 

0.7-2.1 

- 

- 

0.9 ± 0.3 

Progenitor Mass (M©) 

14.0 ± 3.0 

12-13 

8-15.5 

- 

15 ±4.2 

Pre-SN Radius (R©) 

450±112 

230-600 

- 

400-600 

250 ± 70 

AfjVi (Mq) 

0.019 ±0.002 

0.02 ±0.01 

- 

- 

0.013 ± 0.001 

Plateau Duration (Days) 

~ 85 

~ 50 

- 

- 

94 ± 7 

Distance Assumed (Mpc) 

9.57 ± 0.7 

9.6 ±0.7 

9.1 ± 1.0 

9.1 

- 

Distance Measured (Mpc) 

- 

- 

- 

- 

9 0^® ^ 


‘‘Mass quoted is the ZAMS mass of the progenitor, elsewhere it is the final progenitor mass immediately before explosion. 


Bersten et al. ~| |2011[ [Anderson et al. ||2014[ ). 


Applying the i‘‘ model to the early-time data, we esti¬ 
mated the shock breakout epoch of SN 2013ej to be MJD 
56496.9 ± 0.3 days; however, the validity of a model 
in the context of SNe IIP at very early time remains 
to be studied thoroughly. The late time light curve of 
SN 2013ej shows a broken decline behavior in all BVRI 
bands. While M^i = 0.019±0.001 Mq, derived from the 
bolometric lightcurve before the break poi nt at -|-183d. 


is consistent with the prev ious studies (e.g. Bose et al. 
2015at[Huang et al. |2015D, this is perhaps overestimated 
by 50%. Beyond the break point, the slope is shallower, 
close to the expected rate of decay from ^®Co of 0.01 mag 
day“^, resulting in Mmi = 0.013 ± 0.001 M©. The char- 
acterstic time scale of trapping beyond the break point is 
much longer than that found earlier, possibly suggesting 
that SN 2013ej had not completely transitioned to the 
nebular phase before -|-183d, or there was some excess 
flux from CSM interaction with ejecta or other source 
before that time. 

Collecting multiband photometry from U through K 
for a few well-observed SNe from the literature, we estab¬ 
lish a calibration relation between the UBVRI pseudo- 
bolometric flux and the UBVRIJHK bolometric flux, 
which may reach 2% precision. By performing a com¬ 
posite calibration, we showed that ROTSE or KAIT un¬ 
filtered measurements together with B — V information 
may also be sufficient to derive the bolometric luminosity 
with high precision for SN IIP. We also present a pseudo- 
bolometric BVRI calibration using a linear relation for 
unfiltered photometry. Given the position of SN 2013ej 
in relation to more typical SNe IIP and more rapidly de¬ 
clining SNe IIL, it would be interesting to explore the 
consistency of the bolometric calibration described here 
with a broader range of SNe II. 

Even though SN 2013ej is mostly normal spectroscop¬ 
ically, the strong early appearance and subsequent evo¬ 
lution of the Si II A6355 line is rather unusual. Fe II 
might also have appeared somewhat early. The veloc¬ 
ity evolution of weak and strong ions resembles usual 
SN IIP behavior, but the flat H 1 velocity profile is con¬ 
sistently high among the SN IIP population. The corre¬ 
lation of vpe II and uh /3 as observed in a better-sampled 
study is also justified by SN 2013ej. One could expect 
somewhat different behavior with two components (pho- 
tospheric and high velocity) of H I lines, but we do not 
see any HVFs in our spectra. SN 2013ej adds a valuable 
UV spectrum to the early-time UV sample, supporting 
the spectral homogeneity around -1-10 d after explosion. 
It is clear from the available sample that such homogene¬ 
ity is not justified at earlier epochs. This also signifies 
the need of early time data. 


By performing an EPM analysis of SN 2013ej, in com¬ 
bination with SN 2002ap, we estimated the distance to 
the host galaxy M74 to be 9.0to'6 Mpc. Using the cal¬ 
ibrated bolometric flux for unfiltered ROTSE data for 
SN 2013ej, we obtain a distance of 9.Tig;® Mpc, consis¬ 
tent with the previous derivation. Various physical pa¬ 
rameters derived here and the fin ding s from other studies 
of SN 2013ej are listed in Table [^ Generally, the val¬ 
ues derived here from simple approximation models are 
consistent with other findings in the literature. Since 
we took the peak as the advent of plateau, we note that 
the expelled mass is likely overestimated while the radius 
could be underestimated. 
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56535.31 

12.61 

(0.02) 

56536.25 

12.63 

(0.03) 

56537.31 

12.66 

(0.01) 

56537.24 

12.64 

(0.01) 

56538.19 

12.66 

(0.01) 

56539.32 

12.70 

(0.01) 

56539.22 

12.69 

(0.01) 

56540.31 

12.72 

(0.02) 

56540.23 

12.71 

(0.01) 

56541.23 

12.75 

(0.01) 

56542.21 

12.76 

(0.01) 

56543.31 

12.77 

(0.01) 

56543.23 

12.78 

(0.01) 

56549.36 

12.88 

(0.01) 

56549.23 

12.86 

(0.01) 

56563.21 

13.07 

(0.01) 

56563.16 

13.04 

(0.01) 

56565.28 

13.11 

(0.02) 

56565.20 

13.07 

(0.01) 

56567.24 

13.11 

(0.01) 

56568.23 

13.14 

(0.01) 

56568.17 

13.14 

(0.01) 

56569.24 

13.14 

(0.01) 

56569.11 

13.09 

(0.02) 

56570.26 

13.19 

(0.02) 

56570.19 

13.17 

(0.02) 

56571.27 

13.18 

(0.02) 

56571.19 

13.19 

(0.02) 

56572.24 

13.20 

(0.02) 

56572.21 

13.21 

(0.02) 

56573.16 

13.22 

(0.02) 

56574.27 

13.25 

(0.02) 

56575.22 

13.26 

(0.02) 

56575.16 

13.24 

(0.02) 

56576.20 

13.27 

(0.02) 

56576.18 

13.27 

(0.02) 

56577.37 

13.27 

(0.02) 

56578.12 

13.30 

(0.02) 

56586.38 

13.59 

(0.04) 

56586.34 

13.54 

(0.03) 

56587.18 

13.55 

(0.03) 

56587.18 

13.59 

(0.02) 

56588.18 

13.61 

(0.02) 

56588.17 

13.61 

(0.02) 
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56589.19 

56589.15 
56590.24 

56590.16 

56591.19 

56591.15 

56592.13 

56593.17 

56593.15 

56594.27 

56599.16 

56602.19 

56602.27 

56603.27 
56605.33 
56605.31 

56606.17 

56606.21 

56607.22 

56608.12 

56608.13 
56616.09 
56616.09 

56617.13 

56617.20 

56618.17 

56627.27 

56627.13 

56628.22 

56630.21 

56630.23 

56631.16 

56631.23 

56650.20 
56651.19 

56651.16 

56653.21 

56670.16 

56672.15 

56672.13 

56673.13 

56673.15 

56674.12 

56675.13 

56675.18 

56676.14 

56676.16 

56677.11 

56678.12 

56678.11 

56679.12 

56679.17 

56682.21 

56683.17 

56687.15 


13.63 

(0.02) 

13.62 

(0.01) 

13.68 

(0.01) 

13.66 

(0.01) 

13.72 

(0.01) 

13.72 

(0.01) 

13.77 

(0.01) 

13.88 

(0.03) 

13.82 

(0.03) 

13.96 

(0.04) 

14.83 

(0.02) 

15.24 

(0.03) 

15.23 

(0.03) 

15.38 

(0.04) 

15.44 

(0.05) 

15.44 

(0.04) 

15.47 

(0.04) 

15.44 

(0.04) 

15.55 

(0.06) 

15.65 

(0.07) 

15.51 

(0.04) 

15.55 

(0.08) 

15.80 

(0.07) 

15.58 

(0.08) 

15.73 

(0.04) 

15.56 

(0.13) 

15.74 

(0.09) 

15.72 

(0.15) 

15.70 

(0.24) 

15.78 

(0.11) 

15.95 

(0.07) 

15.85 

(0.07) 

16.00 

(0.07) 

16.17 

(0.12) 

16.09 

(0.08) 

16.23 

(0.11) 

16.52 

(0.14) 

17.04 

(0.50) 

16.99 

(0.34) 

16.57 

(0.15) 

16.30 

(0.07) 

16.74 

(0.14) 

16.63 

(0.12) 

16.98 

(0.22) 

16.68 

(0.11) 

16.88 

(0.19) 

17.03 

(0.19) 

16.58 

(0.19) 

16.51 

(0.12) 

16.96 

(0.16) 

16.87 

(0.14) 

16.58 

(0.07) 

17.00 

(0.21) 

16.61 

(0.11) 

16.78 

(0.15) 


TABLE 12 

BVRI PHOTOMETRY OF SN 2013EJ FROM KONKOLY 
Observatory. Magnitudes are in the Vega-system, and 

UNCERTAINTIES ARE GIVEN INSIDE PARENTHESES. 


MJD 

B 

(mag) 

V 

(mag) 

R 

(mag) 

I 

(mag) 

56505.5 

12.64 

(0.05) 

12.57 

(0.02) 

12.49 

(0.02) 

12.43 

(0.02) 

56507.5 

12.61 

(0.01) 

12.52 

(0.01) 

12.39 

(0.01) 

12.36 

(0.01) 

56510.5 

12.66 

(0.07) 

12.49 

(0.07) 

12.35 

(0.02) 

12.29 

(0.01) 

56511.5 

12.75 

(0.08) 

12.51 

(0.04) 

12.31 

(0.01) 

12.26 

(0.01) 

56512.5 

12.78 

(0.10) 

12.52 

(0.03) 

12.31 

(0.01) 

12.24 

(0.02) 

56513.5 

12.78 

(0.06) 

12.48 

(0.04) 

12.27 

(0.04) 

12.21 

(0.01) 

56516.5 

12.99 

(0.08) 

12.53 

(0.02) 

12.31 

(0.01) 

12.20 

(0.02) 

56517.5 

13.05 

(0.04) 

12.54 

(0.02) 

12.30 

(0.02) 

12.19 

(0.01) 

56521.5 

13.38 

(0.06) 

12.61 

(0.02) 

12.32 

(0.01) 

12.17 

(0.01) 

56522.5 

13.39 

(0.05) 

12.64 

(0.01) 

12.34 

(0.01) 

12.19 

(0.02) 

56535.5 

14.15 

(0.06) 

13.06 

(0.01) 

12.61 

(0.03) 

12.42 

(0.04) 

56537.5 

14.22 

(0.01) 

13.11 

(0.01) 

12.65 

(0.01) 

12.43 

(0.01) 

56540.5 

14.42 

(0.07) 

13.19 

(0.07) 

12.70 

(0.01) 

12.50 

(0.01) 

56542.5 

14.46 

(0.03) 

13.23 

(0.01) 

12.74 

(0.01) 

12.51 

(0.01) 

56543.5 

14.49 

(0.02) 

13.26 

(0.01) 

12.76 

(0.01) 

12.52 

(0.01) 

56555.4 

14.89 

(0.12) 

13.48 

(0.01) 

12.95 

(0.06) 

12.66 

(0.02) 

56563.4 

15.01 

(0.02) 

13.59 

(0.03) 

13.03 

(0.01) 

12.77 

(0.01) 

56568.6 

15.09 

(0.01) 

13.69 

(0.02) 

13.11 

(0.01) 

12.84 

(0.01) 

56573.6 

15.29 

(0.10) 

13.77 

(0.01) 

13.19 

(0.02) 

12.91 

(0.01) 

56578.3 

15.34 

(0.06) 

13.88 

(0.01) 

13.27 

(0.01) 

12.98 

(0.01) 

56591.4 

15.93 

(0.04) 

14.37 

(0.01) 

13.67 

(0.06) 

13.40 

(0.01) 

56592.5 

15.94 

(0.05) 

14.41 

(0.01) 

13.73 

(0.03) 

13.43 

(0.01) 

56597.3 

16.65 

(0.03) 

15.10 

(0.01) 

14.28 

(0.02) 

13.99 

(0.01) 

56604.4 

17.65 

(0.04) 

16.24 

(0.01) 

15.25 

(0.02) 

14.92 

(0.02) 

56628.2 

17.95 

(0.07) 

16.68 

(0.01) 

15.63 

(0.02) 

15.36 

(0.01) 

56629.2 

17.96 

(0.06) 

16.70 

(0.01) 

15.68 

(0.03) 

15.43 

(0.01) 
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TABLE 13 

Sloan g’r’i’z’ photometry of SN 2013ej from Baja 
Observatory, Hungary. Magnitudes are in the AB-system, 
AND UNCERTAINTIES ARE GIVEN INSIDE PARENTHESES. 


MJD 

g' 

(mag) 

r' 

(mag) 

i' 

(mag) 

z' 

(mag) 

56504.08 

12.54 

(0.06) 

12.61 

(0.03) 

12.74 

(0.09) 

13.01 

(0.17) 

56505.07 

12.49 

(0.06) 

12.58 

(0.04) 

12.63 

(0.09) 

12.94 

(0.23) 

56510.08 

12.56 

(0.10) 

12.44 

(0.02) 

12.56 

(0.02) 

12.68 

(0.05) 

56512.09 

12.60 

(0.09) 

12.42 

(0.02) 

12.55 

(0.03) 

12.66 

(0.05) 

56513.06 

12.54 

(0.04) 

12.40 

(0.02) 

12.51 

(0.04) 

12.62 

(0.06) 

56516.06 

12.64 

(0.04) 

12.39 

(0.02) 

12.49 

(0.03) 

12.61 

(0.05) 

56520.06 

12.79 

(0.04) 

12.41 

(0.02) 

12.50 

(0.02) 

12.59 

(0.05) 

56521.06 

12.84 

(0.04) 

12.43 

(0.02) 

12.51 

(0.03) 

12.61 

(0.05) 

56522.08 

12.91 

(0.04) 

12.45 

(0.02) 

12.51 

(0.05) 

12.63 

(0.13) 

56535.00 

13.48 

(0.10) 

12.78 

(0.04) 

12.90 

(0.13) 

12.76 

(0.07) 

56536.05 

13.52 

(0.05) 

12.75 

(0.02) 

12.75 

(0.03) 

12.75 

(0.04) 

56539.07 

13.61 

(0.05) 

12.81 

(0.02) 

12.83 

(0.02) 

12.79 

(0.06) 

56539.97 

13.63 

(0.03) 

12.83 

(0.02) 

12.83 

(0.02) 

12.77 

(0.05) 

56541.92 

13.69 

(0.05) 

12.87 

(0.02) 

12.87 

(0.03) 

12.80 

(0.06) 

56543.10 

13.75 

(0.03) 

12.88 

(0.03) 

12.92 

(0.06) 

12.85 

(0.19) 

56543.95 

13.75 

(0.03) 

12.91 

(0.02) 

12.91 

(0.02) 

12.83 

(0.04) 

56552.88 

13.97 

(0.08) 

13.05 

(0.04) 

13.02 

(0.05) 

12.92 

(0.06) 

56558.91 

14.06 

(0.06) 

13.14 

(0.04) 

13.13 

(0.04) 

12.97 

(0.06) 

56559.90 

14.04 

(0.08) 

13.14 

(0.03) 

13.12 

(0.03) 

12.99 

(0.08) 

56568.88 

14.25 

(0.03) 

13.25 

(0.03) 

13.26 

(0.04) 

13.09 

(0.06) 

56575.98 

14.37 

(0.05) 

13.41 

(0.04) 

13.47 

(0.05) 

13.19 

(0.07) 

56576.95 

14.61 

(0.09) 

13.44 

(0.03) 

13.44 

(0.04) 

13.26 

(0.06) 

56578.88 

14.45 

(0.08) 

13.50 

(0.04) 

13.51 

(0.03) 

13.29 

(0.06) 

56586.87 

14.73 

(0.12) 

13.65 

(0.04) 

13.67 

(0.05) 

13.41 

(0.10) 

56590.95 

14.93 

(0.06) 

13.79 

(0.02) 

13.82 

(0.03) 

13.53 

(0.05) 

56592.96 

15.25 

(0.13) 

13.95 

(0.04) 

14.00 

(0.03) 

13.65 

(0.04) 

56594.88 

15.36 

(0.03) 

14.14 

(0.02) 

14.19 

(0.02) 

13.79 

(0.04) 

56596.03 

15.78 

(0.13) 

14.37 

(0.06) 

14.45 

(0.04) 

14.02 

(0.10) 

56598.01 

15.98 

(0.06) 

14.67 

(0.03) 

14.84 

(0.05) 

14.19 

(0.11) 

56603.90 

16.88 

(0.06) 

15.39 

(0.02) 

15.54 

(0.04) 

14.86 

(0.07) 

56612.81 

16.88 

(0.34) 

15.70 

(0.14) 

15.73 

(0.10) 

15.22 

(0.15) 

56626.97 

17.21 

(0.10) 

15.81 

(0.04) 

16.04 

(0.06) 

15.42 

(0.11) 

56627.89 

17.19 

(0.06) 

15.78 

(0.03) 

16.02 

(0.04) 

15.40 

(0.07) 

56628.75 

17.36 

(0.08) 

15.87 

(0.04) 

16.03 

(0.04) 

15.45 

(0.08) 

56629.90 

17.38 

(0.08) 

15.83 

(0.04) 

16.07 

(0.04) 

15.40 

(0.07) 

56636.81 

17.28 

(0.11) 

15.94 

(0.05) 

16.21 

(0.06) 

15.48 

(0.08) 

56650.78 

17.64 

(0.10) 

16.11 

(0.03) 

16.47 

(0.06) 

15.82 

(0.12) 

56660.84 

17.95 

(0.13) 

16.32 

(0.05) 

16.78 

(0.08) 

15.88 

(0.14) 

56663.78 

17.86 

(0.12) 

16.39 

(0.06) 

16.78 

(0.09) 

15.99 

(0.12) 
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TABLE 14 KAIT Unfiltered and KAIT + Nickel BVRI photometry of 
SN 2013ej. Magnitudes are in Vega system and uncertainties are given inside 

parenthesis. 


M,JD 

Open CCD 
(mag) 

B 

(mag) 

V 

(mag) 

R 

(mag) 

I 

(mag) 

56498.45 

13.30 (0.10) 





56499.44 

12.95 (0.04) 





56504.39 

12.29 (0.06) 

12.65 (0.06) 

12.68 (0.03) 

12.53 (0.03) 

12.52 (0.03) 

56507.51 

12.23 (0.02) 

12.62 (0.03) 

12.54 (0.03) 

12.38 (0.03) 

12.37 (0.04) 

56509.52 

12.24 (0.03) 

12.65 (0.10) 

12.57 (0.05) 

12.38 (0.02) 

12.37 (0.03) 

56515.51 

12.21 (0.04) 

12.91 (0.04) 

12.51 (0.04) 

12.26 (0.03) 

12.17 (0.03) 

56516.51 

12.21 (0.02) 





56519.48 

12.29 (0.02) 

13.22 (0.04) 

12.61 (0.03) 

12.32 (0.02) 

12.20 (0.02) 

56520.51 

12.33 (0.01) 

13.33 (0.04) 

12.64 (0.02) 

12.33 (0.02) 

12.23 (0.02) 

56522.50 

12.36 (0.02) 

13.38 (0.04) 

12.64 (0.02) 

12.30 (0.01) 

12.20 (0.02) 

56523.45 

12.41 (0.03) 

13.50 (0.03) 

12.68 (0.02) 

12.37 (0.02) 

12.23 (0.02) 

56524.50 

12.40 (0.03) 

13.54 (0.04) 

12.70 (0.02) 

12.34 (0.02) 

12.21 (0.03) 

56525.34 

12.45 (0.03) 

13.86 (0.09) 

12.87 (0.04) 

12.44 (0.03) 

12.33 (0.02) 

56527.38 

12.52 (0.02) 

13.82 (0.04) 

12.85 (0.03) 

12.47 (0.04) 

12.32 (0.04) 

56530.49 

12.57 (0.03) 

13.91 (0.04) 

12.92 (0.02) 

12.49 (0.02) 

12.32 (0.03) 

56532.53 

12.62 (0.05) 





56533.43 

12.67 (0.03) 

14.11 (0.04) 

13.01 (0.02) 

12.58 (0.02) 

12.40 (0.02) 

56534.49 

12.64 (0.03) 

14.06 (0.03) 

13.00 (0.02) 

12.55 (0.02) 

12.37 (0.02) 

56536.40 

12.72 (0.03) 

14.23 (0.03) 

13.08 (0.02) 

12.67 (0.02) 

12.46 (0.02) 

56537.45 


14.25 (0.12) 

13.10 (0.05) 

12.62 (0.02) 

12.35 (0.21) 

56539.41 

12.77 (0.03) 

14.30 (0.03) 

13.12 (0.02) 

12.65 (0.02) 

12.46 (0.03) 

56540.41 

12.78 (0.02) 

14.33 (0.03) 

13.15 (0.01) 

12.67 (0.02) 

12.48 (0.02) 

56541.36 

12.84 (0.02) 





56542.42 

12.82 (0.03) 

14.41 (0.03) 

13.19 (0.02) 

12.70 (0.02) 

12.48 (0.02) 

56543.38 

12.87 (0.04) 

14.48 (0.04) 

13.22 (0.02) 

12.75 (0.02) 

12.53 (0.03) 

56545.48 

12.86 (0.02) 

14.47 (0.04) 

13.23 (0.02) 

12.74 (0.02) 

12.54 (0.03) 

56546.41 

12.88 (0.03) 

14.54 (0.04) 

13.30 (0.03) 

12.79 (0.04) 

12.56 (0.04) 

56548.39 

12.92 (0.01) 

14.59 (0.03) 

13.31 (0.01) 

12.82 (0.01) 

12.57 (0.02) 

56549.40 

12.92 (0.01) 

14.65 (0.03) 

13.33 (0.01) 

12.84 (0.01) 

12.61 (0.01) 

56550.43 

12.95 (0.03) 

14.60 (0.03) 

13.32 (0.02) 

12.82 (0.02) 

12.59 (0.02) 

56551.45 

12.97 (0.03) 

14.62 (0.03) 

13.33 (0.02) 

12.83 (0.02) 

12.60 (0.02) 

56552.29 

12.98 (0.03) 

14.79 (0.08) 

13.44 (0.04) 

12.91 (0.04) 

12.70 (0.04) 

56553.36 

13.02 (0.03) 





56559.37 

13.08 (0.02) 

14.90 (0.04) 

13.50 (0.03) 

12.93 (0.04) 

12.72 (0.04) 

56562.31 

13.13 (0.01) 

14.96 (0.03) 

13.56 (0.01) 

13.01 (0.01) 

12.77 (0.02) 

56563.40 

13.15 (0.03) 





56564.37 

13.15 (0.04) 

14.93 (0.04) 

13.58 (0.02) 

13.01 (0.02) 

12.76 (0.02) 

56566.49 

13.14 (0.01) 

15.03 (0.05) 

13.59 (0.02) 

13.02 (0.01) 

12.78 (0.02) 

56567.38 

13.18 (0.02) 





56569.29 

13.21 (0.01) 





56570.39 

13.21 (0.02) 

15.06 (0.03) 

13.65 (0.02) 

13.07 (0.02) 

12.82 (0.02) 

56571.31 

13.22 (0.03) 





56572.35 

13.25 (0.02) 

15.14 (0.04) 

13.70 (0.03) 

13.11 (0.03) 

12.86 (0.04) 

56573.35 

13.22 (0.04) 





56574.35 

13.27 (0.01) 

15.18 (0.05) 

13.74 (0.02) 

13.15 (0.01) 

12.89 (0.01) 

56575.37 

13.30 (0.03) 





56576.33 

13.33 (0.04) 





56577.38 

13.34 (0.03) 

15.27 (0.04) 

13.77 (0.03) 

13.16 (0.02) 

12.94 (0.03) 

56579.33 

13.36 (0.02) 





56580.32 

13.41 (0.03) 

15.38 (0.04) 

13.87 (0.02) 

13.25 (0.02) 

13.02 (0.02) 

56581.24 

13.44 (0.02) 





56587.32 

13.61 (0.03) 

15.65 (0.05) 

14.10 (0.02) 

13.45 (0.02) 

13.20 (0.03) 

56588.33 

13.64 (0.03) 





56589.32 

13.70 (0.04) 





56590.32 

13.72 (0.02) 

15.79 (0.05) 

14.23 (0.02) 

13.56 (0.02) 

13.32 (0.02) 

56591.31 

13.74 (0.04) 





56592.30 

13.89 (0.03) 





56596.29 

14.30 (0.02) 

16.49 (0.06) 

14.91 (0.02) 

14.12 (0.01) 

13.85 (0.02) 

56597.27 

14.45 (0.02) 





56598.28 

14.63 (0.03) 





56599.26 

14.83 (0.02) 

17.00 (0.08) 

15.54 (0.03) 

14.62 (0.03) 

14.35 (0.03) 

56600.27 

15.03 (0.02) 





56601.26 

15.16 (0.02) 





56603.31 

15.38 (0.05) 





56604.27 

15.38 (0.03) 





56605.28 

15.42 (0.03) 





56606.27 

15.45 (0.02) 

17.60 (0.09) 

16.36 (0.05) 

15.27 (0.04) 

14.95 (0.04) 

56608.22 

15.59 (0.06) 





56615.27 

15.55 (0.03) 

17.70 (0.29) 

16.55 (0.09) 

15.39 (0.03) 

15.06 (0.04) 

56618.22 

15.63 (0.02) 

17.86 (0.22) 

16.57 (0.05) 

15.45 (0.02) 

15.17 (0.03) 

56619.26 

15.64 (0.02) 





56620.25 

15.65 (0.03) 





56621.28 

15.68 (0.03) 

17.83 (0.18) 

16.65 (0.05) 

15.49 (0.03) 

15.24 (0.04) 
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56622.23 

15.72 (0.04) 





56624.23 

15.72 (0.03) 

18.06 (0.22) 

16.65 (0.06) 

15.49 (0.03) 

15.24 (0.04) 

56625.25 

15.74 (0.03) 





56626.21 

15.72 (0.03) 





56627.23 

15.77 (0.02) 

18.18 (0.25) 

16.70 (0.05) 

15.57 (0.03) 

15.31 (0.03) 

56628.25 

15.78 (0.03) 





56629.25 

15.81 (0.02) 





56630.20 

15.83 (0.02) 

18.13 (0.18) 

16.72 (0.05) 

15.61 (0.03) 

15.37 (0.03) 

56631.18 

15.88 (0.03) 





56632.17 

15.89 (0.04) 





56634.24 

15.90 (0.02) 





56635.15 

15.92 (0.03) 





56636.17 

15.97 (0.04) 





56641.20 

15.98 (0.03) 

18.17 (0.30) 

16.88 (0.08) 

15.78 (0.03) 

15.54 (0.04) 

56643.20 

16.02 (0.06) 





56645.24 

16.03 (0.04) 





56647.23 

16.06 (0.03) 





56648.21 

16.08 (0.03) 

18.06 (0.20) 

17.13 (0.08) 

15.83 (0.03) 

15.63 (0.03) 

56649.20 

16.10 (0.02) 





56650.18 

16.11 (0.04) 





56651.19 

16.15 (0.02) 

18.53 (0.36) 

17.05 (0.08) 

15.91 (0.03) 

15.68 (0.04) 

56653.18 

16.18 (0.03) 





56655.15 

16.23 (0.02) 

18.23 (0.30) 

17.12 (0.11) 

15.98 (0.04) 

15.80 (0.05) 

56656.17 

16.24 (0.04) 





56658.15 

16.30 (0.04) 

18.19 (0.21) 

17.15 (0.07) 

16.01 (0.03) 

15.80 (0.04) 

56660.15 

16.39 (0.06) 





56662.17 

16.32 (0.02) 





56668.12 

16.46 (0.05) 





56669.12 

16.41 (0.04) 





56673.15 

16.52 (0.05) 





56674.11 

16.48 (0.04) 

18.06 (0.25) 

17.47 (0.10) 

16.22 (0.05) 

15.99 (0.05) 

56676.16 

16.55 (0.05) 





56677.13 

16.54 (0.04) 





56679.13 

16.56 (0.04) 

18.53 (0.41) 

17.54 (0.11) 

16.33 (0.05) 

16.11 (0.06) 

56682.12 

16.62 (0.06) 





56684.13 

16.61 (0.04) 

18.63 (0.30) 

17.66 (0.14) 

16.39 (0.05) 

16.21 (0.05) 

56708.12 

16.87 (0.05) 





56710.13 

16.94 (0.05) 





Nickel Photometry 

56505.4258 


12.37 (0.02) 

12.52 (0.01) 

12.43 (0.01) 

12.40 (0.01) 

56507.4883 


12.39 (0.01) 

12.50 (0.01) 

12.37 (0.01) 

12.32 (0.01) 

56523.3945 


13.40 (0.03) 

12.66 (0.01) 

12.36 (0.01) 

12.21 (0.01) 

56527.3398 


13.77 (0.04) 

12.79 (0.01) 

12.45 (0.01) 

12.26 (0.01) 

56531.3125 


13.90 (0.02) 

12.94 (0.01) 

12.50 (0.01) 

12.32 (0.01) 

56535.2969 


14.14 (0.02) 

13.05 (0.01) 

12.60 (0.01) 

12.39 (0.01) 

56539.3242 


14.29 (0.02) 

13.17 (0.01) 

12.69 (0.01) 

12.45 (0.01) 

56541.3008 


14.38 (0.02) 

13.19 (0.01) 

12.72 (0.01) 

12.49 (0.01) 

56545.3008 


14.51 (0.02) 

13.27 (0.01) 

12.79 (0.01) 

12.52 (0.01) 

56548.2695 


14.59 (0.01) 

13.33 (0.01) 

12.84 (0.01) 

12.57 (0.01) 

56552.3516 


14.67 (0.02) 

13.41 (0.01) 

12.89 (0.01) 

12.63 (0.01) 

56555.3203 


14.80 (0.02) 

13.45 (0.01) 

12.92 (0.01) 

12.68 (0.02) 

56559.3086 


14.92 (0.01) 

13.52 (0.01) 

12.98 (0.01) 

12.72 (0.01) 

56562.3555 


15.00 (0.01) 

13.55 (0.01) 

13.00 (0.01) 

12.73 (0.01) 

56569.2930 


15.02 (0.03) 

13.54 (0.04) 

13.13 (0.04) 

12.78 (0.01) 

56573.3047 


15.27 (0.01) 

13.74 (0.01) 

13.17 (0.01) 

12.91 (0.01) 

56575.2656 


15.30 (0.01) 

13.78 (0.01) 

13.18 (0.01) 

12.91 (0.01) 

56578.3086 


15.42 (0.01) 

13.74 (0.01) 

13.25 (0.01) 

12.98 (0.01) 

56581.2812 


15.52 (0.01) 

13.94 (0.01) 

13.31 (0.01) 

13.04 (0.01) 

56585.3867 


15.67 (0.01) 

14.07 (0.01) 

13.43 (0.02) 

13.19 (0.03) 

56589.2656 


15.87 (0.01) 

14.24 (0.01) 

13.57 (0.01) 

13.28 (0.01) 

56599.2617 


17.07 (0.02) 

15.62 (0.01) 

14.70 (0.01) 

14.41 (0.01) 

56601.2734 


17.39 (0.01) 

15.97 (0.01) 

14.99 (0.01) 

14.67 (0.01) 

56605.2695 


17.68 (0.01) 

16.33 (0.01) 

15.28 (0.01) 

14.94 (0.05) 

56609.2656 


17.79 (0.02) 

16.44 (0.01) 

15.36 (0.01) 

15.04 (0.01) 

56626.2031 


18.07 (0.01) 

16.72 (0.01) 

15.60 (0.04) 

15.31 (0.08) 

56631.2148 


18.13 (0.03) 

16.79 (0.01) 

15.68 (0.01) 

15.39 (0.01) 

56635.2266 


18.17 (0.03) 

16.84 (0.01) 

15.73 (0.01) 

15.46 (0.01) 

56642.3359 


18.27 (0.16) 

16.78 (0.06) 

15.89 (0.03) 

15.53 (0.03) 

56645.2578 


18.41 (0.07) 

16.96 (0.02) 

15.89 (0.01) 

15.61 (0.01) 

56665.1875 


18.47 (0.04) 

17.32 (0.03) 

16.21 (0.02) 

15.93 (0.03) 

56670.1992 


18.66 (0.05) 

17.36 (0.02) 

16.27 (0.01) 

16.03 (0.01) 

56675.1914 


18.64 (0.04) 

17.44 (0.03) 

16.27 (0.06) 

16.03 (0.09) 

56678.1719 


18.66 (0.02) 

17.51 (0.02) 

16.32 (0.04) 

16.08 (0.08) 

56689.1250 


18.81 (0.03) 

17.63 (0.02) 

16.52 (0.01) 

16.25 (0.01) 

56705.1211 


18.96 (0.08) 

17.86 (0.03) 

16.71 (0.02) 

16.47 (0.04) 

56709.1250 


18.78 (0.06) 

17.78 (0.04) 

16.74 (0.02) 

16.51 (0.03) 

56724.1367 


19.17 (0.12) 

17.91 (0.04) 

16.88 (0.02) 

16.68 (0.03) 

56849.4336 


19.96 (0.19) 

19.27 (0.08) 

18.56 (0.04) 

18.34 (0.08) 

56852.4062 


20.65 (0.34) 

19.20 (0.08) 

18.46 (0.05) 

18.09 (0.06) 
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56856.4141 


19.74 (0.15) 

19.21 (0.11) 

18.35 (0.06) 

18.41 (0.08) 

56894.3438 


20.31 (0.07) 

19.88 (0.04) 

19.21 (0.03) 

18.49 (0.07) 

56898.4609 


20.27 (0.05) 

20.02 (0.05) 

19.08 (0.05) 

18.73 (0.08) 

56916.4062 


20.49 (0.10) 

20.27 (0.07) 

19.55 (0.04) 

18.78 (0.06) 

56923.4453 


20.55 (0.06) 

20.09 (0.06) 

19.70 (0.05) 

19.03 (0.08) 

56931.3242 


20.86 (0.11) 

20.40 (0.07) 

19.84 (0.05) 

19.61 (0.15) 

56934.3633 


20.77 (0.14) 

20.21 (0.07) 

19.69 (0.07) 

19.21 (0.09) 

56937.2617 


21.06 (0.38) 

20.36 (0.20) 

19.59 (0.10) 

19.02 (0.10) 

56948.2344 


20.56 (0.11) 

20.41 (0.09) 

19.82 (0.06) 

19.70 (0.13) 

56958.2734 


20.72 (0.10) 

20.74 (0.10) 

20.13 (0.07) 

19.54 (0.12) 
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TABLE 15 

Optical and NIR Photometry of SN 2013ej from RATIR. Magnitudes are in AB-System and Photometric uncertainties are 

GIVEN INSIDE PARENTHESIS 


MJD 

r 

i 

Z 

Y 

J 

H 


(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

56500.0 

13.0048 (0.02) 

13.2938 (0.02) 

13.3672 (0.02) 

13.3388 (0.02) 

13.8606 (0.05) 

14.1553 (0.07) 

56502.0 

12.7652 (0.02) 

12.9600 (0.02) 

12.8912 (0.02) 

13.1330 (0.02) 

13.4210 (0.05) 

13.8766 (0.07) 

56503.0 

12.6862 (0.02) 

12.8271 (0.02) 

12.8809 (0.02) 

13.0350 (0.02) 

13.3518 (0.05) 

13.6459 (0.07) 

56504.0 

12.5889 (0.02) 

12.7816 (0.02) 

12.8486 (0.02) 

12.9499 (0.02) 

13.3239 (0.05) 

13.5421 (0.07) 

56505.0 

12.5562 (0.02) 

12.7271 (0.02) 

12.8202 (0.02) 

12.9174 (0.02) 

13.2535 (0.05) 

13.4951 (0.07) 

56506.0 

12.4911 (0.02) 

12.6626 (0.02) 

12.7144 (0.02) 

12.9752 (0.03) 

13.2002 (0.05) 

13.4347 (0.07) 

56508.0 



12.7077 (0.02) 

12.8287 (0.03) 

13.1325 (0.05) 

13.3669 (0.07) 

56509.0 

12.4615 (0.02) 

12.6214 (0.02) 

12.6969 (0.02) 

12.8081 (0.03) 



56510.0 

12.4368 (0.02) 

12.6400 (0.02) 

12.6499 (0.02) 

12.7811 (0.03) 

13.0620 (0.05) 

13.3554 (0.07) 

56511.0 

12.4306 (0.02) 

12.6635 (0.02) 

12.6288 (0.02) 

12.7200 (0.02) 

13.0414 (0.05) 

13.3042 (0.07) 

56512.0 


12.6110 (0.02) 

12.5960 (0.02) 

12.7521 (0.03) 

12.9975 (0.05) 

13.3020 (0.07) 

56517.0 

12.4202 (0.02) 

12.5405 (0.02) 

12.5384 (0.02) 

12.6508 (0.03) 

12.8858 (0.05) 

13.1608 (0.07) 

56521.0 


12.6113 (0.02) 

12.5129 (0.02) 

12.6143 (0.02) 

12.8340 (0.05) 

13.0919 (0.07) 

56523.0 

12.5081 (0.02) 

12.5448 (0.02) 

12.6007 (0.02) 

12.5571 (0.02) 

12.8797 (0.05) 

13.0437 (0.07) 

56541.0 

12.8768 (0.02) 

12.8752 (0.02) 

12.7664 (0.02) 

12.7870 (0.03) 



56552.0 

13.0548 (0.02) 

13.0406 (0.02) 

12.8758 (0.02) 

13.0518 (0.03) 

13.1447 (0.03) 

13.3756 (0.04) 

56559.0 

13.1300 (0.02) 

13.1293 (0.02) 

12.9241 (0.02) 

13.1574 (0.03) 

13.2173 (0.03) 

13.4510 (0.04) 

56571.0 

13.2975 (0.02) 

13.3096 (0.02) 

13.0648 (0.02) 

13.2636 (0.03) 

13.3389 (0.03) 

13.6176 (0.04) 

56578.0 

13.4078 (0.02) 

13.4114 (0.02) 

13.1695 (0.02) 

13.5447 (0.03) 

13.4835 (0.03) 

13.7575 (0.04) 

56585.0 

13.5672 (0.02) 

13.6186 (0.02) 

13.2659 (0.02) 

13.6321 (0.02) 

13.6367 (0.03) 

13.9024 (0.04) 

56588.0 

13.7055 (0.02) 

13.7697 (0.02) 

13.4106 (0.02) 


13.7585 (0.03) 

14.0406 (0.04) 

56589.0 

13.7021 (0.02) 

13.8036 (0.02) 

13.4152 (0.02) 

13.8345 (0.02) 

13.7986 (0.03) 

14.1202 (0.04) 

56593.0 

13.9689 (0.02) 

14.0149 (0.02) 

13.5742 (0.02) 

13.9862 (0.02) 

13.9705 (0.03) 

14.1668 (0.05) 

56595.0 

14.1545 (0.02) 

14.2544 (0.02) 

13.7627 (0.02) 

14.2009 (0.02) 

14.1914 (0.03) 

14.4074 (0.05) 

56600.0 

15.0154 (0.02) 

15.2823 (0.02) 

14.5393 (0.02) 

15.0502 (0.02) 

15.1220 (0.03) 

15.2142 (0.05) 

56604.0 

15.3536 (0.02) 

15.5905 (0.02) 

14.8402 (0.02) 

15.4730 (0.02) 

15.3862 (0.03) 

15.5139 (0.03) 

56608.0 

15.4652 (0.02) 

15.6747 (0.03) 

14.9609 (0.03) 

15.6256 (0.04) 

15.5659 (0.04) 

15.6812 (0.04) 

56623.0 

15.6803 (0.02) 

15.9268 (0.02) 

15.2347 (0.03) 

15.8464 (0.04) 

15.8741 (0.06) 

15.9175 (0.14) 



